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Laminar Boundary-Layer Oscillations and Transition on 

a Flat Plate ' 

By Galen B. Schubauer and Harold K. Skramstad 

This is an account of an investigation conducted at the National Bureau of Standards, 
with the cooperation and financial assistance of the National Advisory Committee for 
Aeronautics, in which oscillations were discovered in the laminar boundary layer along a 
flat plate. These oscillations were found during the course of an experiment in which transi- 
tion from laminar to turbulent flow was being studied on the plate as the turbulence in the 
wind stream was being reduced to unusually low values by means of damping screens. The 
first part of the paper deals with experimental methods and apparatus, measurements of 
turbulence and sound, and studies of transition. A description is then given of the manner 
in which oscillations were discovered and how they were found to be related to transition, 
and then how controlled oscillations were produced and studied in detail. The oscillations 
are shown to be the velocity variations accompanying a wave motion in the boundary 
layer, this wave motion having all of the characteristics predicted by a stability theory 
based on the exponential growth of small disturbances. A review of this theory is given. 
The work is thus experimental confirmation of a mathematical theory of stability which 
had been in the process of development for a period of approximately 40 years, mainly by 
German investigators. 



I. Introduction 

When any fluid flows past solid boundaries, the 
fluid near the surface is retarded by friction. 
This retarded layer is known as the boundary 
layer. While there are many kinds of boundary 
layers, the two main classes are laminar layers and 
turbulent layers. In the vast majority of cases 
of practical interest, a laminar layer will be present 
on the upstream portion of a body while a tur- 
bulent layer will be found on the after part of a 
body. Although a turbulent layer may cover the 
greater part of a body, the laminar layer may be 
regarded as the more natural one, because it must 
exist by virtue of the viscosity of the fluid. A 
turbulent layer is always formed from a laminar 
layer through a break-down of laminar flow re- 

1 This paper was originally published by the National Advisory Committee 
for Aeronautics in April L943asan Advance Confidential Report. The ma- 
terial has recently been released for Unrestricted publication. The introduc- 
tion has been revised somewhat for the benefit of readers not familiar with 
boundarydayer flow. 



suiting from an unstable condition. The nature 

of this instability and the 4 fundamental causes of 
transition to turbulent (low have been the subject 
of many investigations, and theory and experi- 
ment have not been in agreement. The present 
experiment was undertaken to learn more about 
this controversial subject. One of the simplest 
cases was selected for the investigation, namely, 
the boundary layer on a thin flat plate with air 
flowing parallel to the surface. When the pres- 
sure is uniform, the laminar boundary layer on 
such a surface is known as the Blasius type. 

1. Symbols and abbreviations 

#=distance from leading vd^v of Hat plate 
y = distance from surface of flat plate 
Uo= mean velocity outside boundary layer 
£7= mean velocity at a point in boundary layer 
^^instantaneous .r-component of fluctuation 
velocity 
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u 
v' 
w' 



v== instantaneous ?/-component of fluctuation 
velocity 

w= instantaneous fluctuation velocity perpen- 
dicular to Uq and parallel to surface of flat 
plate 

- root-mean-square values of u, v, and w 

c r =w&ve velocity 

fi r =2wf, where /= oscillation frequency 

a~2ir/\, where X= wavelength 
Pi= amplification coefficient 

p= density 

n= viscosity 

^=/x/p= kinematic viscosity 

q= dynamic pressure 

5 = boundary-layer thickness 
$*= boundary-layer displacement thickness 

Vvx 
Yf for Blasius velocity distribu- 
tion 
5* = 0.341 8 relation used by Tollmien and 
Schlichting [2, 4, and 5] 
R= U o*/v= Reynolds number 
Rx= Uox/v=x-^ejno\ds number 

R=l.l2^R x for Blasius velocity dis- 
tribution 

A=-y = Karman-Pohlhausen parameter. 

Past experiments have all indicated that dis- 
turbances in the form of velocity fluctuations, 
either in the oncoming stream or in the boundary 
layer itself, are the direct cause of transition. 
The accumulated evidence favored the view that 
transition depended on the extent to which a 
boundary layer was disturbed; that is, that a 
laminar boundary layer could withstand a certain 
amount of disturbance without becoming un- 
stable, but that it became unstable when the 
disturbances were large enough. It was assumed 
by some that the velocity variations had to be 
large enough to cause intermittent separation of 
the layer from the surface to produce transition. 
Following this point of view, G. I. Taylor [l] 2 ad- 
vanced a theory of transition based on the vary- 
ing pressure gradients accompanying turbulence 
in the surrounding flow outside the boundary 
laver. 



2 Figures in brackets indicate the literature references at the end of this 
paper. 



In spite of this experimental evidence, a sta- 
bility theory based on a quite different concept 
was favored by some mathematicians, mainly 
belonging to the German school. According to 
this theory, small disturbances in velocity of any 
wave length lying within a certain region would 
be amplified in the boundary layer, whereas dis- 
turbances of shorter or longer wave length would 
be damped. Thus stability would depend on 
the wave length rather than the magnitude of a 
disturbance. 

Important contributions to this theory were 
made by Tollmien [2, 31, and Schlichting [4 to 7]. 
There was a little experimental evidence to sup- 
port the theory found in flows known to be of the 
unstable type: namely, flow in a divergent channel 
[8], in a wake at low speeds [9], and in an acous- 
tically sensitive jet [10]. However, this evidence 
was discounted on the basis of inherent instability, 
and stability theory was discredited by those who 
were inclined to trust experiment rather than 
theory. 

Such was the state of the transition problem 
at the outset of the present investigation. The 
original purpose of the investigation was not to 
settle any controversy, but rather to find out 
where transition would occur on a flat surface and 
to learn as much as possible about transition 
when the turbulence of the surrounding steam 
was reduced to levels lower than had hitherto 
been obtained. It was a fortunate choice of 
experiment, for the discovery of laminar boundary- 
layer oscillations, made possible by the smoothness 
of the flow, was the first conclusive evidence of the 
validity of stability theory based on the growth of 
small disturbances. The subsequent study, des- 
cribed herein, of the waves producing the oscilla- 
tions completely verified the theory. 

The failure of previous investigations was due 
in large measure to disturbances that were 
initially large enough to cause transition, these 
disturbances coming in general from the relatively 
high turbulence in the surrounding stream. The 
first known observations of velocity fluctuations 
in a laminar boundary layer were made by Dryden 
[11]. These were regarded as disturbances im- 
pressed on the boundary layer by the turbulence 
in the outside stream. They were irregular and 
produced transition without showing any of the 
characteristics required by stability theory. 
Nikuradse [12] attempted to test the stability 
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theory by producing sinusoidal fluctuations in 
the boundary layer near the leading edge of a 
flat plate in water. His results were inconclusive, 
partly because the artificial fluctuations, as well 
as the general disturbances, were too large, and 
partly because he observed only the effect on 
the location of transition. 

The investigation of transition on a flat plate 
at low turbulence was suggested by Hugh L. 
Dryden, and the authors acknowledge his as- 
sistance and many valuable suggestions in con- 
nection with the experimental program. 

II. Apparatus and Methods 
1. Wind Tunnel 

The present investigation was conducted in the 
4K-foot wind tunnel at the National Bureau of 
Standards. The general layout of the tunnel is 
shown in figure 1. The fiat plate was located 



a honeycomb and resulted in great improvement 
in the steadiness of the stream. Further reduc- 
tion of turbulence was obtained by installing 
damping screens in the settling chamber. With- 
out screens, the turbulence in the test section was 
0.27 percent at 100 feet per second and, with 
seven screens, was 0.032 percent. 

2. Flat Plate 

The flat plate was a commercial sheet of alumi- 
num ){ inch thick, 4% feet wide, and 12 feet long. 
The leading edge, which was symmetrical and 
pointed, was formed by the intersection of circular 
arcs tangent to the original surface 4 inches from 
the leading edge. The surface was left in its 
original condition and had a mirrorlike finish, 
marred slightly by small scratches barely percep- 
tible to the finger tips. Waviness was quite per- 
ceptible when viewed near grazing incidence. A 
surface gage showed variations from approxi- 




Figure 1. — Elevation view of 4 ] /2-ft wind tunnel. 



vertically in the test section with the leading edge 
6 feet from the upstream end. In order to reduce 
vibration, the test section of this tunnel is sup- 
ported directly from the foundation and is joined 
to the rest of the circuit by only a sponge-rubber 
seal at each end. The area reduction from settling 
chamber to test section is 7.1:1. The guide vanes 
ahead of the settling chamber were made finer 
than the others in order to reduce the scale of the 
turbulence and to permit as much reduction in 
turbulence through decay as possible. Since di- 
rectional fluctuations in the horizontal plane were 
found to be large, closely spaced straight en ers 
were placed at right angles to and on top of the 
fine guide vanes. This combination was in effect 



mately 0.01 to 0.02 inch over distances from 1 to 
2 feet. These were bends in the plate evidently 
produced by the rolling process. 

The plate was bolted along the top and the 
bottom edges to one side of 3-inch steel channels, 
which were in turn bolted to the floor and to ths 
ceiling of the tunnel. The flanges away from the 
plate served as rails at the floor and at the ceiling 
for a carriage on which all exploring apparatus 
was mounted. The carriage was simply a steel 
plate % inch thick and (> inches wide, with runners 
and guides at the ends in contact with the rails. 
This plate was thus parallel to the wind and 3 
inches from the flat plate. All explorations in 
the boundary layer were made 25 inches or more 
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ahead of the carriage to avoid the local pressure 
field. The carriage was propelled by hand from 
outside of the tunnel by a sprocket and chain. 

3. Pressure Distribution 

In order to control the pressure gradient along 
the flat plate, the cross section of the working 
chamber was varied by adjustable auxiliary walls 
on the vertical sides of the chamber. These walls 
were aluminum sheets }{ 6 inch thick and 22 inches 
wide extending from opposite the trailing edge of 
the plate to a distance 4 feet ahead of the leading 
edge. The sheets, mounted on screws, spaced 9 
inches, and threaded through the tunnel walls, 
were adjustable to and from the sides of the 
tunnel. In order to prevent discontinuity at the 
top and at the bottom edges, the sheets were 
backed by a stretched rubber diaphragm; the 
combination of aluminum sheet and rubber 
formed a continuous flexible wall that could be 
warped or displaced as desired. The maximum 
range of displacement of these auxiliary walls 
was about 4 inches. 

By this means the pressure could be made to 
rise, to fall, or to remain constant along the plate. 
(See curves A, B, and C, fig. 30.) For most of the 
experiment the walls were set for zero pressure 
gradient. 

The provision for warping the walls into a bulge 
or a hollow did not prove so useful as was antici- 
pated because the effect of the warp usually ex- 
tended over too great a distance to give the de- 
sired result. When it was desired to produce a 
sharp pressure rise or fall, other devices were used. 
For example, an airfoil extending from floor to 
ceiling near the plate produced a sharp pressure 
fall followed by an abrupt rise. A pressure fall 
without a rise was easily produced by proper 
blocking of the stream. 

Although measurements were taken on only 
one side of the flat plate, the auxiliary walls on 
each side of the tunnel were given the same shape 
for symmetry. Some asymmetry nevertheless 
existed because of the presence of the carriage and 
the measuring equipment on only one side of the 
plate. Although this resulted in greater blocking 
of the stream on the working side, it was found 
that some excess blocking here was necessary to 
direct the flow at a slight angle to the leading 
edge of the plate so that the stagnation point was 
displaced slightly from the sharp edge to the 



working side of the plate. The transition point 
was unaffected by the position of the stagnation 
point so long as this point was not displaced to the 
opposite side of the plate. The slight excess 
blocking was therefore maintained at all times to 
prevent directional variations accompanying the 
turbulence in the stream from ever carrying the 
stagnation point to the opposite side. 

4. Measurement of Pressure and Transition 

A small pitot-static head was arranged as shown 
in figure 2 with the impact tube in contact with 
the surface and the static tube parallel to the 




^*3! 
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Figure 2. — Pitot-static-surj'ace tube used to measure pres- 
sure distribution and position of transition. 

stream }i inch from the surface. This combination 
was carried on an arm attached to the carriage 
28 inches to the rear and could slide fore and aft 
with fork-like guides and impact tube always in 
contact with the surface. Both tubes were made 
from thin-wall nickel tubing with 0.04-inch outside 
diameter. The impact tube was flattened on the 
end to form a slit 0.007 inch wide, the center 
of which was 0.006 inch from the surface when 
the tube was in contact. The static tube had four 
0.008-inch holes drilled through the wall, S 
diameters from the closed end. 

Only the static tube of this instrument was used 
when pressure distributions were measured and 
both tubes were used when transition points were 
determined. Traverses were always horizontal 
and usually made along the center line of the 
plate. The static pressure was thus measured 
K inch from the surface, which for most of the 
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surface corresponded to a position outside the 
boundary layer. A large variation of pressure 
normal to the surface was found near the leading 
edge where the pressure distribution was deter- 
mined by the shape of the leading edge. Also, 
pressure effects caused by waviness of the surface 
depended on the distance from the surface. Since 
the former variations were confined to the first 
6 inches from the leading edge and the latter were 
small , variations in pressure in the }£-inch distance 
normal to the surface were neglected. 

The detection of transition by the surface-tube 
method depends on the variation with x of impact 
pressure near the surface, which corresponds 
closely to the variation of shearing stress at the 
surface. This impact pressure decreases with 
distance from the leading edge to the beginning 
of transition, then rises through the transition 
region, and again falls in the turbulent region. In 
the present work, the point of minimum pressure 
was taken as the beginning of transition, and 
the point of maximum pressure following the rise 
was taken as the point where turbulence was 
fully developed. 

5. Boundary-Layer Thickness 

Velocity distributions aeross the boundary layer 
were determined by traversing normal to the 
surface with a flattened impact tube similar to 
the surface tube. This tube was carried on an arm 
similar to the one shown with the hot-wire 
anemometer attached in figure 3. Both the 







Figure 3.- — Apparatus for traversing normal to the surface. 
Uot-wire head attached. 

oiicrometer screw, A, and the fulcrum, B, were 
held in contact with the surface of the plate by 
the rigid arm extending rearward to the carriage, 
and motion of the impact tube to and from the 
surface was obtained by rocking the arm about 
the fulcrum by means of the screw. The velocity 



distributions were used to determine 5*. When 
the pressure gradient was zero, the measured 
values of 5* were in agreement with those cal- 
culated by the Blasius formula 



:i - 72 V 



vx 



Whenever the value of 5* was needed for pi esenting 
results obtained at zero pressure gradient, the 
value calculated by this formula was used. 

Tollmien and Schlichting, [2,4,5] used a constant 
in the above formula equal to 1.73 for computing 
<5*. More recent values of the constant are 1.7207 
given by Dryden [11] and 1.7208 given by Gold- 
stein [13]. A value of 1.72 is therefore used in 
the present work. 

6. Kot-Wire Turbulence Equipment 

A variety of apparatus, such as amplifiers, 

bridges, potentiometers, oscillators, oscillographs, 
and numerous types of hot-wire anemometer, 
comes under the heading of hot-wire equipment 
for measuring turbulence. Amplifiers with attend- 
ant power supplies and hot-wire anemometers are 
usually special equipment, whereas the other 
equipment is of standard commercial design. 
Amplifiers used in turbulence work and circuits to 
compensate for the thermal lag of hot wires have, 
however, been standardized to a certain extent, 
and their essential features are described m the 
literature. Two amplifiers, which were designed 
and built by W. C. Mock, Jr., were used in the 
present investigation. One of the amplifiers is 
described in reference 14; the other was a newer 
and more portable type with about the same 
frequency response but using capacitance com- 
pensation instead of inductance compensation. 
Both types of compensation and the requirements 
to be met by each type are described in reference 
15. The theory of compensation is given in 
reference 16. When properly compensated, the 
over-all response of wire and amplifiers was 
uniform from 3 to about 2,000 cycles per second. 
In all cases hot-wire anemometers must be 
calibrated by measuring the voltage 4 across the 
wire at known air speeds; or, if the anemometer is 
sensitive to direction changes, it must be cali- 
brated by measuring the voltage at several 
angles to the wind. For this purpose a poten- 
tiometer is necessary. A Wheatstone bridge 
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must be used to obtain the resistance of the wire 
at air temperature. Throughout the present 
work, both in calibration and in use, the heating- 
current was held constant and only the tempera- 
ture and the voltage across the wire were allowed 
to vary. This condition was maintained during 
rapid velocity fluctuations by a sufficiently large 
choke coil in series with the wire. Amplifiers 
were calibrated by applying to them a known 
alternating voltage from an oscillator. By use of 
the calibration data, velocity fluctuations could 
be calculated from the amplified and properly 
compensated voltage fluctuations across the hot 
wire. A cathode-ray oscillograph was used for 
visual observation of fluctuations and also for 
making photographic records with a moving-film 
camera. When observations were simultaneously 
taken on two wires, an electronic switch was used 
in conjunction with the oscillograph. 

7. Hot-Wire Anemometers 

The hot-wire anemometers used here may be 
divided into two classes: (1) Those used in the 
free stream to measure turbulence, (2) those used 
in the boundary layer to study oscillations, tran- 
sition phenomena, and turbulence. Anemometers 
of class 1 were designed for high sensitivity and 
for freedom from vibrational effects. Anemom- 
eters of class 2 were designed for working near a 
surface with a minimum of interference and with 
as little vibrational motion relative to the surface 
as possible. Both types used platinum wire 
0.00025 inch in diameter, obtained by etching the 
silver coating from Wollaston wire. The etched 
wire was soft-soldered to prongs made either from 
fine sewing needles or from fine copper wire. 
Platinum wires of this size were adequately sensi- 
tive in lengths no greater than Ke inch and, in 
the range of operating conditions, had time 
constants ranging from 0.0005 to 0.002 second. 

Anemometers of class 1 were made with a 
single wire normal to the wind for measuring u f 
or with two wires set at an angle to the wind in 
the form of an X for measuring v' and w'. Heads 
of these two types are shown in figure 4. Both 
have platinum wires (not visible in the photo- 
graph) across the tips of the prongs; the four- 
prong head has the two wires forming an X. In 
the X -arrangement the wires lie as nearly as 
possible in a plane without touching and subtend 
an angle of 60 degrees, with each wire at an angle 





Figure 4. — Hot-wire heads for the measurement of free- 
stream turbulence. 

of about 30 degrees to the wind. Braces of silk 
thread cemented across the prongs about }i inch 
from the tips were necessary to prevent wire 
breakage. All wires were about % 6 inch long. 
These arrangements of wires take advantage of 
the directional characteristics of a wire in order 
to obtain sensitivity to the component desired. 
The sensitivity of a wire to direction is known to 
depend on the angle to the wind, being zero when 
normal and again when parallel. Sensitivity to 
u is a maximum when the wire is normal. Since 
v and w are small, their principal effect is to cause 
fluctuation in direction of the stream when added 
vectorially to U and the effect on the magnitude 
of the instantaneous velocity is insignificant. A 
wire normal to the wind responds thus to u- 
fluctuations rather than to v or w. 

An X -arrangement of two identical wires, each 
making the same angle to the wind, with voltage 
opposed produces a resultant voltage change only 
when the wires are differentially cooled. This 
arrangement is then insensitive to u, but highly 
sensitive to direction changes in the plane of the 
wires and hence to v or to w depending on the 
orientation of the plane, When fluctuations are 
large, isolation of single components in this way 
is not complete and voltage changes are not 
exactly proportional to velocity changes. If the 
fluctuations do not exceed 5 percent of the mean 
speed, errors from these sources are believed to 
be less than 1 percent. 

The mountings for the hot-wire heads are shown 
in figure 5. The central member to which the 
head is attached is held within the cylindrical tube 
by a rubberband suspension and in this way is 
isolated from vibration from the tunnel walls. 
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Vibrational effects were not troublesome after a 
sufficiently "soft" suspension of this sort was used. 
The mounting for the X-wires contains an angle- 
changing device for calibrating the wires. 



fei,to>_:.:.u.l_ 




Figure 5.— Rubber-suspension-type mountings with hot- 
wire heads attached. 

In the boundary layer only the ^/-component 
was measured, except when one attempt was made 
to discover whether the boundary-layer oscilla- 

I tions had a w-eomponent. A study of v was con- 
sidered desirable but no appropriate hot-wire head 
of sufficiently small dimensions was available for 
the purpose. For anemometers of class 2, there- 
fore, most attention was given to hot-wire heads 
with the wire normal to the wind and parallel to 

I the surface. The wire length was usually 0.04 
inch. For traversing along the surface with the 
wire at a fixed distance, the wire was attached at 
the tips of prongs extending about 1/4 inch forward 
from a small celluloid sled, which was held against 
the surface by wire springs from a brass tube at- 
tached to the carriage. A sled assembly with two 
pairs of prongs for supporting two wires at differ- 
ent distances from the surface is shown in figure 6. 
Another sled was also used with two hot wires 




Figure 6. — Sled-like hot-wire head arranged for two hot 

wires. 
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making a V in a plane parallel to the surface. The 
latter arrangement was sensitive to w. The sled 
was believed to be the best possible type of hot- 
wire head to prevent relative motion between the 
wire and the surface. Interference effects were 
found negligible. 

For traversing normal to the surface, the hot 
wire was supported on the tips of steel sewing 
needles extending into the boundary layer from 
the arm of the traversing apparatus shown in 
figure 3. This arrangement was made sufficiently 
rigid to be free from vibration troubles but was 
not appropriate for traversing parallel to the sur- 
face. A combination of this apparatus with a 
sled was used when it was desired to keep one wire 
at a fixed distance from the surface while moving 
the other in and out from the surface. 

III. Reduction of Stream Turbulence and 
Effect on Transition, Zero Pressure 
Gradient 

In the present investigation, effects of high 
stream turbulence were of little concern. Em- 
phasis was placed rather on reducing the turbu- 
lence as far as practicable in order to study the 
boundary layer when it was as little disturbed 
from the stream as possible. The turbulence was 
successively reduced by placing various numbers 
and combinations of damping screens in the set- 
tling chamber. For each addition of a screen and 
for every recombination of screens, measurements 
were made of u 1 '/U , v'/U , and w' '/U in the work- 
ing chamber and the values of B x at the beginning 
and at the end of transition were determined. 

No correction for wire length according to the 
method given in reference 17 was made, since the 
appropriate scale for v f and w f could not be de- 
termined. If the probable order of magnitude of 
the scale is considered, the correction is small and 
is believed to be unimportant in view of the 
accuracy attainable in the measurement of low 
turbulence. 

Screens of fine mesh and small wire were used 
in order that the turbulence produced by the 
screen itself would decay rapidly. Each screen 
completely covered the cross section of the set- 
tling chamber and, when more than one screen 
was used, a spacing of 6 inches or more was 
allowed between them. Preliminary surveys with 
the hot-wire anemometer showed that the turbu- 
lence was uniform over the usable cross section of 
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the stream and that little or no decrease in turbu- 
lence occurred in the length of the working cham- 
ber. The turbulence could therefore be specified 
by hot-wire measurements made at one point in 
the working chamber. The measurements showed 
that v' and w' were nearly equal to each other but 
generally greater than u f . This noniso tropic con- 
dition is known to be caused by the contraction 
of the stream from the large area of the settling 
chamber to the smaller area of the working cham- 
ber. The more the turbulence was reduced by 
the addition of screens, the less this difference 
became and, finally, near the lowest turbulence 
u' was actually larger than v f and w' . In all cases 




W' 2 )/Uo 



looV'/sd^+v 

Figure 7. — Effect of turbulence on x-Reynolds number of transition. 
Flat plate, zero pressure gradient. 



the turbulence increased with the wind speed, an 
effect that could be partly accounted for by the 
decreasing resistance coefficient of the screens with 
increasing speed. 

As the number of screens was increased, a reduc- 
tion in turbulence was obtained with each addi- 
tional screen but by successively decreasing 
amounts; and the reduction by this means ap- 
peared to have reached a practicable limit when 
six screens were installed. The measured values 
with six screens are given in table 1. It was 
believed at first that a level had been reached 
nearly equal to the turbulence produced by the 
screens themselves. A finer screen than any of 
the others was therefore added downstream from 
the six screens already present. As no screens 
with finer wires than those already used were com- 
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30 
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.0334 


90 
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100 


.037 
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.041 


.0401 


110 


.040 


.044 


.045 


.0430 


120 


. 045 


.046 


.047 


.0459 



mercially available, a screen made from silk bolt- 
ing-cloth was installed. Some father reduction 
was obtained but not so much as was expected. 
No attempt was made at further reduction. The 
measured values of the turbulence at this lowest 
level are given in table 2. The x-Reynolds num- 
bers of transition found during the various stages 
of reduction are shown in figure 7. The highest 
turbulence shown in the figure (0.342 percent) 
was obtained by placing a square-mesh grid of 
/2-inch rope, spaced 6 inches, in the settling cham- 
ber a few inches downstream from the set of six 
wire screens. The transition point found by Hall 
and Hislop [18] for their lowest turbulence is indi- 



cated in figure 7. By definition, their transition 
point corresponds to the beginning of transition 
of the present paper. 

Table 1. — Turbulence with 6 damping screens in settling 

chamber 

Measurements at single station in working chamber 
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Two questions are raised by t! 
suits: (1) Why was the addition 
cloth screen so ineffective in furtl 
turbulence? and (2) Why was ti 
flat plate unaffected by reductio 
below about 0.08 percent? The 
answered in the following sect in 
question is more difficult, and an 
made to answer it only after 
oscillations are considered. (See 



he roregomg re- 
of the bolting- 
ler reducing the 
ansition on the 
n of turbulence 
first question is 
n. The second 
attempt will be 
boundary-layer 
section VIII.) 



IV. Residual Turbulence and Noise 

While the measurements of the lowest turbu- 
lence were in progress, it was noted by the trace 
on the screen of a cathode-ray oscillograph that the 
velocity fluctuations did not have the random 
character usually ascribed to turbulence. Vibra- 
tion of the hot-wire mounting was suspected but 
was ruled out since changes in the softness of the 
mounting failed to change the appearance of the 
trace. Consideration was next given to the effect 
of wind-tunnel noise and led to the startling 
conclusion that a reasonably loud sound could 
involve particle velocities sufficient to produce an 
apparent turbulence of the order of that actually 
found. For example, a calculation based on a 
plane wave with an intensity of 105 decibels above 
a base level of 10~ 9 erg per cm 2 per second showed 
a root-mean-square particle velocity of 0.028 foot 
per second. At a wind velocity of 100 feet per 
second, this is equivalent to an apparent turbu- 
lence of 0.028 percent, while the measured turbu- 
lence was 0.0315 percent. A few measurements of 
intensity with a sound meter connected to a crystal 
microphone showed that the noise level wa"s proba- 
bly between 105 and 110 decibels at 100 feet per 
second. As 110 decibels is equivalent to an appar- 
ent turbulence of 0.051 percent, the noise could 
account for the entire hot-wire reading at this 
speed. At the lower speeds, however, the noise 
intensity appeared insufficient to account for the 
measured turbulence. As the wave pattern in the 
tunnel was very complex, calculations based on 
plane waves can be at best only rough approxima- 
tions. Furthermore, to avoid the direct effect 
of the wind on the microphone, it was necessary 
to measure the sound level at a distance from the 
hot wire. 

An analysis was made by a wave analyzer of the 
spectrum of the noise from the tunnel picked up 



in the tunnel control room by the crystal micro- 
phone. A similar analysis was then made of the 
output from the hot wire in the wind tunnel. In 
both cases the wind velocity was 80 feet per second. 
A comparison of the two spectrums is shown in 
figure 8. As the distribution of intensity with 
frequency is compared here, uncertainties in abso- 
lute intensity are of minor importance. The 
similarity of the curves indicates that, at a wind 
velocity of 80 feet per second, a large part of the 
hot-wire output is derived from noise. 

The values of turbulence in tables 1 and 2 
probably contain a significant contribution from 
noise, especially at the higher velocities. A large 
part of the increase with velocity is probably due 
to noise. Also, a u' greater than v' and w f is 
evidence that not all the contributing fluctuations 
were carried with the stream through (he entrance 
(•one. As most of the noise conies from the pro- 
peller, it is obvious that damping screens can 
reduce only tin 1 true turbulence. On the other 
hand, the addition of screens must increase the 
noise because of the higher propeller speed needed 
to maintain a given wind velocity" through the 
added resistance in the screen. These facts afford 
a satisfactory explanation of the effective limit 
of damping screens in the present tunnel. 

Table 2. — Turbulence with 6 damping screens followed by 
1 bolting cloth in settling chamber 

Measurements al single station in working chamber 



To 


u' 

Co 


v' 

Co 


w' 

Co 




yi/3(« ,2 +t;'2+w'2) 
Co 


fps 


Percent 


Percent 


Percent 


Percent 


30 


0.018 


0.012 


0.009 


0.0135 


40 


.019 


.012 


.010 


. 0142 


50 


.021 


.014 


.015 


. 0109 


60 


. 020 


.018 


.010 


.0204 


70 


.029 


. 020 


.017 


. 0226 


80 


.033 


. 022 


.018 


. 0251 


90 


.037 


. 026 


.020 


.0286 


100 


.010 


.0:9 


.023 


.0315 


110 


' .044 


.034 


. 025 


. 0352 



Figure 8 also suggests some information con- 
cerning the answer to the second question posed 
in section III regarding the effect of turbulence 
on transition. As the spectrum of this apparent 
turbulence is unlike that of isotropic turbulence 
[191, such turbulence will probably affect transi- 
tion in a different manner from isotropic turbu- 
lence. Velocity fluctuations from sound involve 
pressure changes not related to those of turbu- 
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lence and involve motions correlated over great 
distances. As the proportion of noise to turbu- 
lence was increased with each reduction of 
turbulence, anomalous transition effects are to be 
expected. 

V. Laminar Boundary-Layer Oscillations, 
Zero Pressure Gradient 

1. First Evidence of Oscillations 

When the stream turbulence had been reduced 
to nearly its lowest level with the six wire damp- 
ing screens, it was decided to investigate the 
velocity fluctuations in the laminar boundary 
layer itself. For this purpose one of the sled-like 
heads with a wire sensitive to u was arranged to 
slide along the surface with the wire 0.023 inch 
from the surface. With greatly reduced stream 
turbulence, it was expected that the slow irregular 
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50 



- 20 



10 



Records of the oscillations were made by photo- 
graphing the screen of the cathode-ray oscillo- 
graph with a moving film camera. Two sets of 
such records are reproduced as figures 9 and 10. 
Figure 9 shows the progressive development of 
the oscillations. The traces at 6.25 and 6.5 feet 
already show the bursts of large-amplitude oscil- 
lations interspersed with highly disturbed turbu- 
lent motion. The final record at 8 feet indicates 
a completely turbulent boundary layer. The 
frequency of the oscillations may be judged by 
the Ko-second timing dots. The same phenom- 
enon at a lower wind speed and at a greater dis- 
tance from the leading edge of the plate is shown 
in figure 10. Here the amplifier gain (magnifi- 
cation) was adjusted to keep the recorded ampli- 
tude 4 the same from position to position in order 
to bring out the weak as well as the strong oscil- 
lations. Attention is called to the lower oscilla- 
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fluctuations, reported for a much more turbulent 
stream in reference 11, would be correspondingly 
reduced. Such fluctuations were, in fact, almost 
nonexistent; but, as the wire was moved down- 
stream through the boundary layer, a regular 
oscillation appeared, weak at first but with in- 
creasing amplitude as the distance downstream 
increased. Just ahead of the transition, bursts of 
very large amplitude occurred and, at the initial 
point of transition, these bursts were accom- 
panied by a breaking into an irregular high- 
frequency fluctuation characteristic of turbulence. 



tion frequency corresponding in this case to the 
lower speed and to the thicker boundary layer. 

It was thought at first that these oscillations 
might possibly be due to vibration of the wire 
that gave rise to relative motion between the wire 
and the surface. This was soon ruled out since 
changes in the hot-wire head and support failed 
to cause any change in the oscillations. Further- 
more, the oscillations appeared to be definitely 
connected with transition since the zone in which 
they occurred always preceded transition and 
moved with it fore and aft along the plate as the 
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Figure 9. — Oscillograms of u- fluctuations showing laminar 
boundary-layer oscillations in boundary layer of flat 
plate. 

Distance from surface=0.023 in. Uo=80 ft/sec. Time interval between 

dots =1/30 sec. 

wind speed was varied. It could not be said 
with certainty that the oscillations were the cause 
of transition, since there existed the possibility 
that the boundary layer became shock-excited by 
transition occurring a short distance downstream, 
giving rise to an oscillation that was possibly the 
result of transition rather than the cause. The lat- 
ter explanation was ruled out by removing transition 
with an abrupt pressure drop and yet leaving the 
oscillations totally unchanged at an upstream 
position. The alternative conclusion, with no 
evidence to the contrary, was that transition 
resulted from the growing of the oscillations to the 
point where the boundary layer was so highly 
disturbed that transition occurred. 

It was soon found that the hot wire could be 
placed at any cross-section of the bounclaiy layer 
to pick up the oscillations, although the amplitude 
quite evidently varied with position. The fre- 
quency remained the same throughout the section 
except very near the surface where there was a 
suggestion of frequency doubling on the low-veloc- 
ity part of the cycle. This effect at the 4-foot 
position is shown in figure 11. As the hot wire 
responds to negative velocities as though they were 



positive, the doubling was believed to indicate that 
the amplitude of the oscillation was sufficiently 
high to reverse the direction of flow during half of 
each cycle. 

A V-wire, sensitive tow, revealed a w-componenl 
in the oscillations. Little work was done with t his 
wire, and the relative amplitudes of u and w were 
not compared. With the one exception cited here, 
all work was dore with wires sensitive to u. Be- 
cause of experimental difficulties no attempt w r as 
made to detect the ^-component, although it is 
show 7 !) by theoretical considerations in section V, 2, 
that a ^-component must exist. The presence of 
the ^-component thus indicates that the oscilla- 
tions were three-dimensional. 

Because the effects of noise were so much in 
evidence, it might be supposed that the oscillations 
were a resonant acoustic phenomenon arising from 
some frequency in the sound spectrum. Among the 
evidence against this supposition is that presented 
in figure 12, which shows the effect of a tenfold 
increase in turbulence produced by the rope grid 
in the settling chamber downstream from the six 
damping screens. With the higher turbulence, it 
is seen that oscillations are present where they 
were absent before. It thus appears that dis- 



K 8.0 sA/ 



aAaaa 



> 5 ./vA 



z9.0 



vAA^vVW^AA/ww^ 3 1 







Figure 10. — Oscillograms of u-fl actuations showing lamirar 
boundary layer oscillations in boundary layer of flat 
plate. 

Distance from surface=0.023 in. £7o=53 ft/sec. 'rime interval between 
dots = 1/30 sec. 
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Figure 11. — Oscillograms of u-fluctuations showing evidence of frequency doubling near surface. 
Ub= 105 ft/sec. Time interval between dots=l/30 sec. 
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Figure 12. — Oscillograms of u-fluctuations showing effect of stream turbulence on boundary -layer oscillations. 
Distance from surface =0.004 in. Time interval between dots=l/30 sec. 
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Figure 13. — Zone of amplification enclosed by neutral curve, according to Schlichting. 

Theoretical total amplification from branch I to branch II as indicated. Experimental points from frequency of 
oscillations found on oscillograms. Branches I and II of neutral curve meet R-ax\s at /?=«>. 



turbances in the stream, possibly acoustic as well 
as turbulent, give rise to oscillations which are not 
themselves sound waves. 

It was believed almost from the start that the 
oscillations were amplified disturbances, their 
purity resulting from selective amplification of a 
single frequency or, at most, a narrow bund of 
frequencies. Theories had been advanced to 
account for just such an amplification but were 
often ignored in experiment because no conclusive 
evidence in their support had been found. In 
order to test whether these were in reality the 
amplified oscillations predicted by theory, fre- 
quencies were determined from numerous oscillo- 
grams, taken at as many positions and speeds as 
possible and plotted in an appropriate manner on 
a diagram derived from the theoretical results of 
Schlichting. The diagram with the experi- 
mentally determined points is shown in figure 13. 
This figure will be made clear by the theoretical 
discussion in section V, 2, and will be described in 
detail in section V, 3. For the present, it is suffi- 
cient to observe that the experimental points are 
distributed along branch II of the theoretical 
curve, which is the region in which the points 
should fall if the initial disturbance, out of which 
oscillations grow, have received approximately 
their maximum amplification. The oscillations are 
therefore identified with those produced by insta- 
bility of the boundary layer. 

2. Theory of Boundary-Layer Oscillations 

Evidence has been found of oscillations of 
rather high purity resulting mostly from random 



disturbances. The oscillations observed al a 
fixed point suggest the presence of a traveling 
wave in the boundary layer. A theoretical answer 
is sought for the question: How can random dis- 
turbances produce a wave that looks closely 
sinusoidal and has a single predominant frequency? 
The superposit ion of small disturbance velocities 
on a uniform velocity in one direction presents to 
a stationary observer the appearance of wave 
motion. The motion may be sinusoidal or irregu- 
lar. For simplification, all the 4 disturbance 
velocity components are assumed to be sinusoidal 
of the same frequency. in a uniform stream, 
this motion has the attributes of a traveling wave 
with a wave velocity equal to the mean velocity 
of the stream. The amplitude is the disturbance 
velocity with components u, v, and w, and the 
pressure involved is the dynamic pressure of the 
disturbance velocity. This kind of wave is in- 
volved in the application of stability theory to the 
boundary layer. Only two-dimensional motion has 
so far been successfully treated. In a nonuniform 
stream, such as a boundary layer, the wave 
velocily is not known at once but may reasonably 
be assumed (and can also be proved theoretically) 
to be less than the maximum velocity of the stream. 
The characteristics of the wave will depend on the 
characteristics of the stream, such as the velocity 
profile of a boundary layer, and on conditions 
imposed on the mean velocity and on the disturb- 
ance velocities at the boundaries. The theory 
treats the influence of all these conditions and how 
they determine the wavelength, the wave velocity, 
and the damping or the amplification of the wave. 
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In general, the disturbances responsible for 
wave motion in a boundary layer are random in 
character and the resulting wave motion is ir- 
regular. One of the important characteristics 
of the wave is its ability to maintain itself against 
the damping action of viscosity or even to grow by 
absorbing energy from the basic flow. Since this 
characteristic depends on the frequency, certain 
component frequencies of an irregular wave will 
be damped and others amplified as the wave travels 
downstream. It will be seen in section V, 3 
that this process accounts for the relatively high 
purity of the observed boundary-layer oscillations. 

Energy considerations alone form the basis of an 
elementary treatment first used by Osborne 
Reynolds. According to the Reynolds concept, 
growth or decay of a disturbance depends on 
whether energy is transferred to the disturbance 
by absorption of energy from the basic flow or is 
extracted from the disturbance by the damping 
action of viscosity. The flow of energy then 
determines the stability of the boundary layer 
under the action of disturbances. As the kind of 
disturbance must be specified, this method does 
not have general applicability. The flow of energy 
is, however, a fundamental physical concept in- 
volved in all theories even though not explicity 
stated. As an introduction to the more advanced 
theories the energy relation will therefore be 
derived here. The development follows essentially 
that given by Prandtl [20]. 

(a) Basic Equations and Energy Relations 

In general, only two-dimensional flow with two- 
dimensional disturbances has been treated success- 
fully in stability theories. The basic flow is 
assumed to be steady and a function of y only. 
The basic flow in the boundary layer is then 

U=f(y) 
F=0. 

The disturbances are a function of time as well as . 
of x and y and are expressed by 



and products of disturbance velocities and of their 
differential coefficients are neglected, 



d 2 U d 2 u &u 
by 2 bx 2 dy : 



(1) 
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7 p\bx bxj 






r 



U £ =V \^ 2+ ^ 2 )--p(ty + &i,)' (2) 



where p is the pressure produced by the disturb- 
ances, and P is the pressure due to the basic flow. 
The equation of continuity is 



dx dy~ 



(3) 



By subtracting from eq 1 and 2, respectively, the 
corresponding Navier-Stokes equations for u= 
v=p=o, the following two equations are obtained 
in terms of only the disturbances : 



dt 



„, bU /i> 2 u,d 2 u\ 1 dp fA . 

^dy= v (w + wr-pzx (4) 



r C>U 

dx 



dv u dx~ v \dx 2 ^W) p*y () 

Equations 4 and 5 are the fundamental hydro- 
dynamic equations for small disturbances on 
which are based all stability theories herein 
discussed. 

In order to derive an expression for energy 
balance, eq 4 is multiplied by u and eq 5 by v, and 
the resulting equations are added to give the single 
disturbance equation 



d+^)m= 



"["(■ 



d 2 u d 2 u\ . 
dx 2+ dy 2 )^ 



puv 



bU 
dy 






dx 2 ^dy 2 Jj' 



(6) 



u=f 1 (x,y, t) 

v=Mx,y, 0- 

The components of total velocity are therefore 
U+u and v. As u and v are assumed to be small, 
the Navier-Stokes equations become, after squares 



It will be noted that the left-hand member of this 
equation gives the time rate of change of kinetic 
energy, due to the disturbances, of a particle 
moving with the basic flow. Each term on the 
right-hand side of eq 6 is now integrated over the 
region containing the disturbance, as explained in 
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detail in reference [20]. The resulting expression 
is 



<1K 
dt 



~ p f ! uvd ^ dxdy ~"i S(^~^) dxdyf (7) 



where dE/dt now represents the time rate of change 
of kinetic energy, due to the disturbance, of the 
fluid within the region of integration. The first 
term on the right-hand side of eq 7 involves the 
shearing stress puv and represents the rate at 
which energy is absorbed from the basic flow. 
This energy is added to that already present in the 
disturbance. The second term on the right-hand 
side, which is obviously always negative, is the 
rate of dissipation of the energy of the disturbance 
by viscosity. The predominance of one effect 
over the other will determine the net change in the 
energy of the disturbance. 

The first term on the right-hand side of eq 7 
has an interesting physical interpretation. First, 
if u and v differ in phase by 90 degrees, as in the 
usual wave motion, uv will vanish over one 
complete cycle, and no energy can be absorbed 
from the basic flow. Second, as dU/dy is positive, 
uv must be negative, if the disturbance is to 
receive energy from the basic flow. These and 
other conditions are treated in the more advanced 
theory to be considered next. 

(b) Disturbance Equation and Solutions 

Lord Rayleigh [21] appears to have been the 
first to outline a more general mathematical 
theory to which numerous investigators have con- 
tributed. The mathematical difficulties are so 
great that the basic velocity profile (U =/(?/)) was 
usually approximated by a straight line or seg- 
ments of straight lines. This imposed serious 
limitations on the generality of the results. Tiet- 
jens [22] was the first to apply Rayleigh's theory 
to profiles intended to reproduce flow along a wall, 
but his solution suffered from the serious limita- 
tions just mentioned. Important advances have 
been made by Tollmien [2, 3] and Schlichting 
[4 to 7], who successfully applied the theory to 
curved profiles intended to represent actual velo- 
city distributions. Their most complete solutions 
apply to ao approximation of the Blasius distribu- 
tion. Schlichting's work [4, 5] is particularly com- 
plete and his results are used here for comparison 
with experimental results. The following outline 



of the theory is taken largely from references 
[4 and 5] with occasional help from reference [2]. 
The physical basis for the theory, the steps in the 
solution, and the more important results are 
given here. The mathematical aspects of the 
problem are discussed in the original references. 
Even in the original references details of the com- 
putations are lacking and checking of results is 
difficult. 

From the basic eq 4 and 5, the pressure terms 
are eliminated by differentiating eq 4 and 5 with 
respect to y and x, respectively, and subtracting 
the second from the first. The result is a linear 
homogeneous equation in u and v: 
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(8) 



The disturbance velocities may be expressed in 
terms of a stream function \p such that 



u=b\(/ 
dy 



dx 



(9) 



Any periodic disturbance may be represented by a 
Fourier series. As eq 8 is linear and homogeneous 
in u and v, its behavior may be investigated by 
using a single term of the series. This amounts 
to assuming a periodic disturbance whose stream 
function has the form 

f=F(y) exp [i(ax—pt)] = F(y) exp [ia(x-ct)] 9 (10) 

where F (y) represents the initial amplitude of the 
stream function which depends only on y, a=2w/\, 
where X is the wavelength, and t is the time. 
Since £ and hence c are generally complex quan- 
tities, eq 10 may be written 

f=F(y) exp [i(az-(p r +iPt)t)] = 

F(y)exp[ia(x-(c r +ic t )t)], (11) 

where r , the real part of 0, is the angular velocity 
or 2tt/ (/= frequency) ; f , the imaginary part of 0, 
is the coefficient of amplification or damping, 
depending on whether it is positive or negative; 
and c r =p r /a is the velocity with which the phase 
progresses and will be termed simply "wave 
velocity. " 

When eq 8 is written in terms of the stream 
function and substitution is made for \f/ in accord- 
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ance with the second of eq 10, the result is 

(U-c) (F" -a 2 F)-U"F=^(F"" -2ocF" +a i F), 

(12) 

where the primes represent differentiation with 
respect to y. Equation 12 is a homogeneous linear 
differential equation of the fourth order with a 
general solution of the form 



F=C 1 F l + C 2 F 2 +C 3 F,+C i F 4 , 



(13) 



where F u F 2 , F z , and F 4 are the particular solu- 
tions and C u C 2 , C 3 , and C A are the constants of 
integration. The form of these particular solu- 
tions depends on the basic velocity assumed, that 
is, on the boundary -layer profile. 

A solution of eq 12 was first obtained by Toll- 
mien [2] for a velocity distribution intended to 
approximate closely a Blasius distribution. He 
considered the case in which amplification and 
damping were absent, that is, for real values of 
/3 and c, and obtained the so-called "neutral" 
oscillations. Schlichting [4] repeated Tollmien's 
calculation and in addition treated the case for 
small amplification. In a later paper [5] Schlicht- 
ing determined the distribution of amplitude of 
neutral oscillations across the boundary layer and 
investigated the energy balance. Schlichting 's 
work is therefore more complete and more suitable 
for coniDarison with experiment. 

Both Tollmien and Schlichting approximated a 
Blasius distribution by a straight line and a para- 
bola as follows: 
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= 1. 



This gave a sufficiently close approximation to 
the Blasius distribution of U but was not a suffi- 



ciently good approximation for TJ f and U" . The 
latter quantities were therefore taken directly 
from the Blasius solution; only the following terms 
of the series, which obviously apply to small values 
of y/8, were used 



U' 

I 



^1.68[l-3.65(|) 3 ] 



and 



The displacement thickness 

5*= \ S (l-jAdy = 0.M75. 

Schlichting gives 0.3415 for the integral, and this 
value was used to compute the values of <5 in 
figure 21. 

It was assumed that U=f(y) and that V=0. 
This is equivalent to assuming a nonthickening 
boundary layer, that is, constant 8, and is the basis 
for most of the adverse criticism of the theory. 
In a real boundary layer, 5 depends on x, U= 
f(x,y), and W0, but the dependence of U on x is 
small compared with the dependence on y, and 
V is small compared with U. 

As the purpose in the present paper is merely to 
indicate the method of solving the equation, the 
solutions will be indicated in general form. The 
calculations have been carried through to nu- 
merical results for only a few velocity profiles. 

Equation 12 is not readily solved in its entirety, 
and certain approximations must be made that 
involve the friction terms on the right-hand side 
of the equation. The simplest approximation is 
to neglect the friction terms entirely. In order to 
justify this, eq 12 is put in a dimensionless form 
that has the Reynolds number as one of the 
parameters. When the new independent variable 
y/8* is introduced and F=U 8*<t>, the result is 

5SB ( *"*~ 2 " !S *'*" + "' S "* ) ' (14) 

where the primes now denote differentiation with 
respect to y/8* and R is the Reynolds number. 
When R is sufficiently large, the right-hand side 
of equation 14 may be small enough to neglect. 
The effect of the friction terms is actually neg- 
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ligible for all values of ?//<$* except those near the 
surface and near a plane, 3 which is also near the 
surface, where the wave velocity is equal to the 
local stream velocity. It has been shown by 
Tollmien [2] that c must always be less than U ; 
a point will therefore always be found in the 
boundary layer where U=c. Exclusive of this 
point and the surface, two particular solutions, 
fa and fa, may be obtained by solving the friction- 
less equation 
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(15) 



The solutions are expressed as power series de- 
veloped about y n where y c is the position of 
the critical point, U=c. These solutions are 
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(17) 



where TJ" and U c are the values of TJ" and TJ' 
at the critical point. While fa is regular through- 
out the whole layer, fa possesses a singularity at 
the critical point and a correction to eq 17 must 
be found by another approximate solution that 
takes into account the friction at the critical 
point. This second approximation to eq 14, 
valid in the neighborhood of the critical point, 
results from neglecting all but the largest of the 
friction terms. 

For this purpose a new independent variable rj 
is introduced, defined by 



y-i/c 



i^vi) 



1/3 



= «?• 



(18) 



From this it is seen that rj is a dimensionless vari- 
able that is equal to zero at the critical point 
y=y c , positive for y^>y c , and negative for y^iy^ 
Equation 14 is now written with r\ as the inde- 
pendent variable for 0. As this equation will be 
used in the vicinity of y^yc, (U—c)/U can be 
approximated by (U c /U )tri and U" by I '', 
where the primes again denote differentiation with 
respect to y/5*. The equation becomes 



3 As the theory applies only to a boundary layer of constant thickness, t \w 
problem is essentially one-dimensional, and this plane is commonly referred 
to as the critical "point." 



ty' " + <t>"(r,-2it 2 a z 5* 



Je jf, +r,e 2 a 2 5* 2 - i e V5* 4 ) = 0, 



(19) 



where </> is differentiated with respect to 77. As 
e is small, the terms containing e 2 and t A arc 
omitted, and eq 19 becomes 






(19a) 



Equation 19a contains no singularity at U=c, 
and a solution of this equation at the critical point 
indicates the form of fa to be used when y—y^O. 
This form is 

«2=§ii(logj^j+iir)+ 

l-bp-^ + bfi-y^+... (17a) 

By neglecting; the term in eq 19a containing 
t, the homogeneous equation 



v>" " + r,<t>" = 



(20) 



is obtained. The solution of eq 20 contains the 
two particular solutions fa s and fa that are ex- 
pressed in terms of Hankel functions as follows: 

fi P*> (D.(2)ro "I 

* 3 ,4= dn v U2 H ~W 2 \dv, (21) 

J 00 J co ]/3 |_0 _J 

where // (1,a ' ) are Hankel functions of the first 
and second kinds. Equation 20 is valid all the 
way to the surface. It is found that fa diminishes 
very rapidly from the surface outward and is there- 
fore important only near the surface. On the 
other hand, fa increases rapidly with ?/, and will 
not satisfy the requirement that </>=0 at y= «> ; 
hence, C 4 in eq 13 is set equal to zero. The gen- 
eral solution of eq 14 may then be written 



= Cifa + C 2 fa -f Czfa. 



(c) Characteristic-Value Problem 



(22) 



The stability investigation is a characteristic- 
value problem — that is, one in which the boundary 
conditions afford sufficient equations between the 
solutions at the boundaries to determine the values 
of the parameters for which eq 14 is satisfied, 
those values being the characteristic values. At 
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y=8 where U=U = constant, the solutions of eq 
14 are of the simple form 



exp 



exp 



(--*) 

(+«sf.) 



The solution with the positive exponent must be 
ignored as it is infinite at y= <*> . As the outer 
boundary condition, then, <j>' '/<j>=— a8*. Since $ 3 
has already disappeared at the outer boundary, 
the result is 



or 



C l (<t>n+a8*<t> l8 ) + C 2 ($' 28 +a8*cl> 28 )=0 



(23) 



where the subscript 8 indicates the solution at the 
outer boundary and 

$25 = 025 + ^*025 

At the surface, y=0 and <£=<// = 0. The surface 
boundary conditions are then simply 



Ci^Iq + &2 02 Q + Cs<l>3 — 
Ci01 ft +O202 n +C3^3 ft = O, 



(24) 
(25) 



where the subscript denotes solution on the 
surface. 

If <j> in eq 22 does not vanish identically, the 
following determinant of the boundary solutions, 
formed fromeq 23, 24, and 25, must equal zero: 



$15 



$25 

*2« 





<t>L 



Solving the determinant gives 



03 o __ 02 o $l 5 — <Al $2 5 
03 () 02 O $15 — <^1 $25 



(26) 



The left-hand side of eq 26 is a function o 77 ,f 
where 770 is the value of rj at the surface. Since, 
by eq 18, 77 is a function of the parameters ad* 
and R, the left-hand side of eq 26 is also a func- 
tion of these parameters. The right-hand side of 
eq 26 is a function of the parameters ad* and 
c/Uq. Denoting the left-hand side by G and the 
right-hand side by E gives 



O{7 h ) = E( a S*,e/U ). 



(27) 



Equation 27 is complex and may be separated 
into two real equations by separating real and 
imaginary parts. Between these two equations 
any one of the three unknown parameters, a8*, 
R, and c/U , may be eliminated and a relation 
between the other two found. This has been 
done by both Tollmien and Schlichting for c real, 
that is, for oscillations neither amplified nor 
damped. For this case, eq 27 may be written 



O n (VOn) = En (a»5 *, C T / Uq) 



(28) 



where the subscript n denotes neutral oscillations. 
If c r is eliminated, a relation between a8* and R 
is found. This gives the theoretical neutral curve 
shown in figures 20, 25, and 27. By eliminating 
a8*, the relation between c r /U and R, shown in 
figures 26 and 33, is found. Recombination of 
these quantities results in 8 r v/U Q 2 , plotted in 
figures 13, 16, 17, 19, and 24. 

The curves just mentioned extend only to the 
maximum values of R that are of practical interest. 
Tollmien [2] points out that both branches I and 
II of these curves approach the i?-axis as R 
approaches infinity. Branches I and II, there- 
fore, meet at infinity and form a closed curve. 
Amplification occurs within the region enclosed 
by the curve. Damping occurs in all other 
regions. 

The values of the parameters are given in 
table 3. The computations were made by a com- 
bination of analytical and numerical methods. 
It may be pointed out that the values obtained 
by Tollmien and Schlichting show some differ- 
ences. These are probably because of differences 
in the details of the computations. Values 
plotted in the foregoing figures and appearing in 
table 3 are Schlichting's results. 
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Table 3. — Theoretical wave parameters for neutral ocilla- 
tions, according to Schlichting [4] 



Cr 


Branch I 


Branch II 


R 


a8* 


•^X10» 


R 


ad* 


^X10* 


0.200 


7200 


0.077 


2.14 


37,600 


0.149 


0.79 


.250 


3010 


.101 


8.39 


12,000 


.188 


3.92 


.300 


1530 


.129 


25.3 


4,640 


.223 


14.4 


.325 


1150 


.143 


40.4 


3,290 


.238 


23.5 


.350 


893 


.159 


62.3 


2,070 


.251 


42.4 


.375 


736 


.181 


92.2 


1,420 


.264 


69.7 


.400 


633 


.205 


129.5 


1,020 


.274 


107.4 


.420 


606 


.239 


165.7 


713 


.273 


161.0 



Schlichting extended the calculations to include 
small amplification, applying only in the neighbor- 
hood of the neutral curve. The calculations were 
carried out, not by a direct solution of eq 27, but 
by a series development starting from the neutral 
eq 28. Thus 6 (r? ) is found by 



G(Vo) = G(rj 0n ) + (vo—Von) Q' (Von) + 
and E (ad*, c/£7 ) is found by 

E(a5*, clU,)^E n (a n b* } c r /U ) + 



(29) 



(30) 



With the values of G and E given by eq 29 and 
30, amplifications were determined near and within 
the space enclosed by the neutral curve Pi=0. 
Values beyond the range permitted by this method 
were determined by interpolation, by assuming a 
cubic equation for p t at constant R that meets 
both branches of the neutral curve with the known 
slope dftilda at these points. Schlichting's dia- 
gram showing the theoretical amplification is 
given in figures 27 and 28. It will be noted 
that p % has a maximum near the center of the 
zone and falls off toward both branches of the 
neutral curve. In accordance with eq 11, ampli- 
fication of a disturbance is expressed by 



A, 

a; 



rt 2 
-exp Pidt, 



(31) 



where Ai is the amplitude at t h and A 2 is the am- 
plitude at t 2 . With t x taken on branch I of the 
neutral curve and t 2 on branch II, the total ampli- 
fication at constant p r vlU 2 was calculated by 
Schlichting for several parts of the amplification 
zone. Schlichting's values are shown in figures 



13, 19, and 20. The maximum total amplifica- 
tion up to a given value of x, when disturbances 
of all frequencies are present, is determined by 
the equation 



A(?) 

a; 



~-exp( max I p t dt V (32) 



Schlichting found that the relation between this 
maximum total amplification and the ^-Reynolds 
number up to an R x of about 10 6 could be repre- 
sented by the interpolation formula 

^^=0.55 exp (0.555 X10" 5 i? x ). (33) 

The solution thus far has determined several 
characteristics of the wave. For any given value 
of R and the frequency, it is known whether the 
wave will be amplified, damped, or neutral. For 
a given frequency the wave velocity is related to 
U in a known way and is constant over the cross 
section of the layer. In like manner the wave 
length is known in relation to 5*. 

(d) Amplitude, Correlation, and Energy Balance 

The characteristic- value problem gives no in- 
formation about the relative magnitudes of u 
and v nor how they are distributed through the 
boundary layer. How the boundary layer "vi- 
brates," whether as a whole or in parts, remains 
unknown. This is because the boundary condi- 
tions were used to determine only the parameters 
of the problem and not the constants of inte- 
gration. 

In a second treatment of the theory applied to 
the flat plate [5], Schlichting completed the in- 
vestigation by determining the constants in equa- 
tion 22 for neutral oscillations, that is, for /3 t =0. 
He was thus able to calculate relative values of u 
and v and their distribution across the boundary 
layer. The physically important questions of the 
correlation between u and v and the energy balance 
were also answered. The same approximation to 
the Blasius distribution as used earlier was used 
here. 

As absolute values of u and v, which obviously 
depend on the intensity of the impressed disturb- 
ance, are of no interest in such a treatment, one of 
the constants was left undetermined and for con- 
venience was made unity. Thus, 
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and, from the boundary conditions expressed by 
eq 23 to 25, the following values for C 2 and C z are 
obtained: 



$25 



^"^(Si ^o-^hM^ ^o-^o)j 



(34) 



When the values of the boundary solutions and 
their derivatives are known, C 2 and C z may bs 
determined by eq 34. The general solution <£, 
is thereby determined except for an arbitrary 
intensity factor. From the nature of the particu- 
lar solution, <t> has real and imaginary parts, <j> T and 
4>t, which must be known with their derivatives in 
order to determine u and v. 

The stream function ^ in eq 9 and the first 
of eq 10 are used to obtain the real parts of u/U 
and tf/Uo, since F— l7 5*<£, as follows: 



(35) 



yf=K[<l>r COS (aX—P r t) — (j)i sin (ax— fl r t)] 

V 

jj=Ka[<j> r sin (ax—P T t)-\-<t>i COS (ax—/3 r t)] 

where K is an arbitrary intensity factor. The 
corresponding root-mean-square values of u and v 
are 

v! K 






(36) 



(37) 



An expression for the kinetic energy is readily 
obtained by squaring eq 36 and 37 and adding. 
The correlation coefficient defined by 



K= 



uv 



becomes 



K= 



<l>r<t>i—<l>r<l>% 



V(*; 2 +*; 2 )(#+*f) 



(38) 



Values of u'/U , v'/U , and K were calculated 
by Schlichting for neutral oscillations correspond- 
ing to the following two points, one of each branch 
of the neutral curve: 



Branch 


R 


VrV 

TT 2 


ad* 


a 


I 


894 
2070 


62.3X106 
42.4 


0. 159 
. 251 


0. 35 


II 


. 35 



The value of K was chosen to give an average 
value of u'/U =0.05 from to 6. Table 4 gives 
the values of u'/U , v'/Uoj and K calculated 
by Schlichting. 



Table 4. — Theoretical neutral oscillations taken from 
Schliching [5] 

For average u'/Uo=0.05 from y/5=0 to yj8=l 



10 u, 



< 



Branch I; R, 894; «5*, 10.159 



0. 050 
.090 
.130 
.170 
.209 
.250 
.290 
.370 
.451 
.531 
.612 
.693 
.774 
.854 
.935 

1.015 

1.1 

1.2 

1.3 

1.4 

1.5 



0.445 


0. 0546 


.647 


.158 


.804 


.287 


.922 


.434 


1.005 


.608 


1.016 


.795 


1.003 


.973 


0.854 


1.322 


.689 


1.605 


.524 


1.820 


.371 


1.984 


.244 


2.098 


.121 


2.170 


.007 


2.193 


.103 


2.180 


.211 


2.108 


.203 


2.027 


.193 


1.935 


.185 


1.849 


.176 


1.763 


.168 


1.681 



-0. 032 
-.061 
-.155 
-.169 
-.140 
-.109 
-.064 
-.022 
-.0075 



























Branch II; R, 2070; a8*, 0.251 



0.029 
.054 
.074 
.105 
.157 
.209 
.250 
.290 
.370 
.451 
.531 
.612 
.693 
.774 
.854 
.935 

1.015 

1.1 

1.2 

1.3 

1.4 

1.5 



0.661 


0.0655 


.937 


.214 


1.038 


.359 


1.108 


.595 


1.127 


1.028 


1.086 


1.445 


0.998 


1.755 


.906 


2.04 


.686 


2.50 


.504 


2.87 


.340 


3.10 


.205 


3.26 


.087 


3.35 


.022 


3.37 


.123 


3.33 


.216 


3.22 


.305 


3.11 


.292 


2.925 


.272 


2.72 


.252 


2.53 


.234 


2.35 


.218 


2.18 



-0. 121 
-.170 
-.191 
-.146 
-.076 
-.010 
+.015 
+.019 





























l 6*=o.3415. 
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It will be seen in table 4 that u'/U =0 at the 
surface and again at a point near y/d= 0.854 for 
branch I and near ?//S = 0.774 for branch II, after 
which it again increases for greater y/8. The 
values of u'/U are plotted in figure 21, where the 
theory is compared with experiment. At the 
points within the boundary layer where u' '/U Q 
vanishes, a 180° phase shift occurs in u, meaning 
that the instantaneous ^-component of the 
oscillation velocity is decreasing on one side of this 
point while it is increasing on the other. No such 
phase shift occurs in v. Table 4 also shows that 
correlation does exist between u and v and that 
uv is negative. This means that energy is drawn 
from the basic flow; and, in accordance with 
eq 7, this energy must be balanced by the dissipa- 
tion if the oscillation is to be neither amplified 
nor damped. Schlichting evaluated both terms in 
eq 7 by integrating from x—0 to x=\ and from 
y=0 to y= co . The net flow of energy was found 
to be very nearly zero, which indicated that the 
oscillations were nearly neutral. This result was 
interpreted as proof of the correctness of the 
solution to the characteristic-value problem. 

A physical explanation of the correlation 
between u and v is attempted by Prandtl [20]. It 
appears that u and v would differ in phase by 90° 
and uv would vanish in consequence if there were 
no phase shift in v produced by the viscosity near 
the critical point. The viscosity is thus responsible 
for the flow of energy to as well as from the 
disturbance. 

(e) Solutions for Other Conditions 

Tollmien [3] has attempted to generalize the 
solution sufficiently to include boundary-layer 
profiles with an inflection point, such as occur with 
an adverse pressure gradient. For such profiles 
the sign of U" changes from positive to negative in 
passing outward from the surface and, at the 
inflection point, U" = 0. The velocity of the 
neutral wave is shown to be equal to the stream 
velocity at the inflection point, and an important 
consequence of this is the absence of a critical point 
at U=c r in the solution of the frictionless eq 15. 
The conclusion is that profiles with an inflection 
point are inherently unstable, although neutral 
oscillations can still exist. These considerations 
suggest that an adverse pressure gradient is 
favorable to amplification. 

Gortler [23] has investigated the effect of wall 



curvature on the stability of the boundary layer 
and finds that velocity profiles are unstable if a 
change of sign occurs in V+U'/r, where / is the 
radius of curvature of the wall. If an inflection 
point (l7" = 0) exists, the point where U"+ ( r '/r=0 
is farther from the surface than the point where 
U" = for a convex wall and nearer to the 4 surface 
for a concave wall. If no inflection point exists, a 
point where U" '+£7' '/r = can occur only on 
a convex wall. It is concluded that convex 
curvature has a destabilizing effect and thai 
concave curvature has a stabilizing effect. 

Solutions for profiles with inflection points are 
given in references 8 to 10. Roscnbrook [8] 
investigated the flow through a diverging channel 
as the channel was towed through still water. 
Streamlines were made visible by imc particles of 
aluminum in the water, and wavelike motions 
accompanied by the formation of eddies could be 
observed as well as the velocity distribution. A 
solution, following Tollmien [3], was obtained for 
S-shape velocity profiles approximating those 
observed. Fair agreement was obtained between 
the observed and the computed wavelengths. 
Hollingdale [9] obtained the solution appropriate 
to the laminar wake formed in one ease 4 behind a 
flat plate and in another behind an airfoil, both 
towed at low speeds through still water. Oscillat- 
ing wakes were observed at certain Reynolds 
numbers and the wavelength and the wave 
velocity were measured. With the appropriate 
interpretation of theoretical results, fair agree- 
ment between calculated and observed wave 
length and wave velocity was obtained. Savic 
[10] obtained an exact solution of the frictionless 
eq 15 for the two-dimensional jet and obtained 
expressions for the neutral wave length and wave 
velocity in terms of characteristics of the jet. 
The theoretical results were compared with 
observations of wave length and wave velocity 
made on an acoustically sensitive jet and good 
agreement was obtained. 

Schlichting [6] has determined the effect of 
centrifugal force by calculating the curves of 
neutral stability for a boundary layer developed 
on the inner wall of a rotating cylinder. Schlicht- 
ing [7] has determined the effect of density 
gradient on the stability of the boundary layer, 
using the approximation to the Blasius distribu- 
tion given in part b of this section. Neutral 
curves outlining the zone of amplification are 
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given. Stability is principally affected by a 
gravitational field in connection with the density 
gradient through the boundary layer. 

The methods of Tollmien and Schlichting have 
been applied by Pretsch [24] to the boundary 
layer on bodies of revolution. The approximate 
forms of the general disturbance equations from 
which the four particular solutions are obtained 
are found to be the same as those for the two- 
dimensional case, the solutions are therefore the 
same. The solutions are applied to the known 
velocity profile at the stagnation point of a sphere 
and the curve of neutral stability is computed. 

3. Naturally Excited Oscillations 

Naturally excited oscillations are those origi- 
nating from the amplification of small disturb- 
ances naturally present in the stream. These 
oscillations were treated in Section V, and their 
frequencies were experimentally obtained and 
plotted in figure 13. The concept of waves run- 
ning downstream through the boundary layer has 
been introduced; the behavior of the oscillations 
in relation to figure 13 can now be discussed. 

One important feature of figure 13 is that the 
values of amplification coefficient are positive 
within the zone enclosed by the neutral curve and 
have values reaching a maximum near the center 
and falling to zero on each boundary. Outside 
the zone f} t is negative. All disturbances, whether 
damped or amplified, travel downstream and 
therefore in the direction of increasing R, since 6* 
increases with x. For example, a fixed speed Z7 
and a disturbance of frequency such that p r v/U 2 
in figure 13 equals 40X10 -6 are considered. 
Such a disturbance, if it began at the leading edge 
of the plate, would be damped as it passed along 
the boundary layer until it reached a point where 
R attained the value 1,180. From this point to 
J?=2,160 the disturbance would be amplified and, 
as it passed beyond this point, would again be 
damped. If A 1 and A 2 are the amplitudes of the 
disturbance at the beginning and at the end of 
the amplification zone, the ratio A 2 /A 1 = 742. In 
other words, the amplitude has increased 742 
times while crossing the zone. This is the total 
amplification given by eq. 31. A disturbance of 
any given frequency will then be most highly 
amplified when it reaches branch II of the neutral 
curve. 



Figure 8 shows that, in experiment, many fre- 
quencies are represented in the initial disturbance. 
The fact that the observed boundary-layer oscil- 
lations contain one predominant frequency must 
be explained. In order to determine whether 
theory and experiment are in agreement, the place 
on the theoretical diagram where the predominant 
frequency should lie must be found. Instead of 
traversing the diagram of figure 13 with a wave 
of a single frequency as in the foregoing illustra- 
tion, waves of many frequencies, covering a wide 
range of values of p r v/U 2 , are supposed to enter 
the amplification zone through branch I and un- 
dergo amplification at a rate determined by their 
frequency as they proceed to the value of R at 
which the observation is made. Initial disturb- 
ances are, of course, present in the boundary layer 
at all values of R and disturbances enter the 
diagram at all points. Those entering this dia- 
gram ahead of branch I are reduced by damping 
below the prevailing level of the initial disturb- 
ances while those entering behind branch I cannot 
receive full amplification. Disturbances occurring 
at points corresponding to branch I have the 
largest initial value and also receive the highest 
amplification and are therefore assumed to be the 
ones whose development should mainly determine 
the final result at the point of observation. It 
will be noted that amplification increases without 
limit as the frequency decreases if no limit is 
placed on the Reynolds number. The experi- 
mental observation, however, is made at a par- 
ticular value of R and hence of R x ; under this con- 
dition it is obvious that one particular value of 
firv/Uo 2 will become more highly amplified than 
any of the others, this one having the maximum 
total amplification given by eq 33 in terms of R x 
and the initial amplitude A x . This frequency 
obviously is the predominant one that should be 
found on the oscillograph records. The remaining 
question now concerns the position of this most 
prominent frequency on the diagram, that is, 
where the experimental points should lie. The 
problem then becomes one of determining along 
what line starting from branch I and running 
parallel to the i?-axis the total amplification is a 
maximum up to a given value of R. By taking 
into account the values of /3 f given by Schlichting, 
it was found that for uniform energy distribution 
with frequency, this line introduced into the dia- 
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gram along branch I intersects the prescribed value 
of R somewhere between the center of the zone 
and branch II. For more accurate results, both 
the initial energy distribution and the damping 
prior to reaching boundary I must be considered. 
The experimental points should lie inside the zone 
near branch II. The agreement with theory is 
not perfect, since many points lie on branch II 
and outside the zone, but is as good as can be 
expected from observations of this sort. 

The process described here is believed to be 
one of considerable importance since it is the proc- 
ess by which laminar-boundary-layer oscillations 
develop from initially random disturbances. In 
order to develop a pure oscillation, initial disturb- 
ances must be so small that a large amount of 
amplification can occur before transition destroys 
the laminar layer. 

It has been pointed out that the oscillations 
appear to be three-dimensional. This is to be 
expected because of the three-dimensional nature 
of the initial disturbances. The theory, however, 
concerns two-dimensional oscillations, and how 
this should affect the agreement with experimenl 
is not known. 

An apparent discrepancy between figures 7 and 
13 must be pointed out. It will be observed in 
figure 13 that oscillations are indicated for R as 
high as 3,120. This corresponds to an R x of 
3.28 X10 6 and, according to figure 7, is in the 
transition region. The oscillograph records actu- 
ally showed intermittent turbulence, but it was 
still possible to read oscillation frequencies from 
non turbulent parts of the records. The limits of 
the transition region are defined statistically by 
the surface-tube method. The discrepancy is 
therefore only apparent. 

In connection with the surface-tube method, 
another effect was noted, namely, that a rise in the 
surface-tube reading was noted a little ahead of 
the point where turbulence was indicated by the 
hot wire. This is believed to be caused by a 
slightly modified velocity distribution arising from 
the shearing stress given by uv in the oscillations. 
The lower boundary of the transition region in 
figure 7 is therefore a little too low, but the dis- 
placement is small and has no significant effect on 
the shape of the curve. 



VI. Artificial Excitation of Boundary-Layer 
Oscillations, Zero Pressure Gradient 

It was soon realized that a study of boundary- 
layer oscillations could be carried out to better 
advantage if they were not caused by accidental 
disturbances occurring in the wind tunnel but 
were produced by a controlled disturbance of 
known amplitude and frequency at some chosen 
position. This was finally accomplished by the 
vibrating ribbon discussed in section VI, 3 . While 
a search was being made for devices of various 
kinds that would accomplish the desired result, 
effects of sound were tried. Interference from 
naturally excited oscillations was minimized by 
making the artificial disturbances much greater 
than the natural disturbances. The bolting-cloth 
damping screen was also installed with the six 
w4re screens, but the addition of this screen had 
no detectable effect on the natural oscillations. 

1. Sound From Loudspeaker 

The effect of sound intensity on the boundary 
layer is of interest because sound is an unavoidable 
source of disturbance both in wind tunnels and 
in free flight of power-driven aircraft. A 25-watt 
loudspeaker was accordingly installed in the top 
of the tunnel at the leading edge of the plate and 
Uh\ by a variable-frequency oscillator through an 
amplifier. The intensity and the frequency were 
controllable, 

Interesting effects were easily demonstrated. 
For example, boundary-layer oscillations could 
be induced at will merely by choosing the right 
frequency for a particular position and speed. 
Transition could be moved 1 or 2 feet ahead of 
its natural position by the right combination of 
intensity and frequency. In general, a random 
noise from the loudspeaker produced similar re- 
sults, but the effect on the oscillations was not so 
marked. These were casual observations with no 
attempt at quantitative measurements. 

When quantitative work was attempted, it be- 
came apparent that the complicated sound field 
in the tunnel was a decided draw-back. Figures 
14 and 15 show the distribution of sound intensity 
along the center of the working chamber 3 inches 
from the plate. Figure 14 shows a standing wave 
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at 48 cycles per second with either a node or a 
loop at the center of the plate depending on the 
type of microphone — that is, whether crystal 
microphone, responding to sound pressure, or 
ribbon microphone, responding to particle velocity. 
This turned out to be the fundamental frequency 
of the 12-foot tube formed by the plate and the 
tunnel walls. The same standing wave could be 
set up without the loudspeaker when the exciting 



frequency came from the eight-blade wind-tunnel 
propeller running at a speed of 360 rpm. Figure 
15 shows the distribution at 60 cycles per second. 
The ribbon microphone and the hot-wire ane- 
mometer show nearly the same distribution, as 
would be expected. These are only a few examples 
of the sound field and the change in distribution 
with frequency that serve to complicate any 
studies with controlled sound intensity. 
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Figure 16.- — Boundary layer oscillations excited by loudspeaker. 

Points show sound frequency giving maximum amplification. Curve outlines theoretical amplification zone. 

(See fig. 13.) 



When the loudspeaker was used to excite oscilla- 
tions, two procedures were tried: (1) The sound 
frequencies that produced oscillations of maxi- 
mum amplitude at chosen velocities and positions 
were determined and (2) a sound frequency and 
position were chosen and the velocity was then 
varied until the oscillations were a maximum. 
By use of the sound frequencies read on the oscil- 
lator together with the measured U and x, 
Prv/Ul and R were calculated and plotted in 
figure 16. The solid line is again Schliehting's 
curve of neutral stability enclosing the zone of 
amplification. Since both procedures amounted 
to finding the condition For maximum total 

200 



amplification, the points should lie within the 
zone near branch II to be in agreement with 
theory. In view of the difficulties encountered, 
the large scattering of the points is not surprising. 
Another application of the Loudspeaker con- 
sisted in directly connecting the amplified output 
of the hot wire to the loudspeaker' so that regen- 
eration could occur. Disturbances picked up by 
the hotwire were thereby i'vd back as sound, which 
in turn reinforced the initial disturbance when 
the proper phase relation existed, causing the 
system to oscillate. The resulting frequency was 
read from a variable-frequency oscillator con- 
nected to the oscillograph to form Lissajou figures 
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Figure 17. — Regeneration frequencies resulting from feedback from hot wire in boundary 

layer to loudspeaker. 

Curve outlines theoretical amplification zone. (See fig. 13.) 
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with the hot-wire signal. The frequency of re- 
generation was determined by the frequency 
response of the wire, the electrical circuit, the 
loudspeaker, and the tunnel as a whole, as well 
as by the preferred boundary-layer frequency. 
Many frequencies could be found with no wind 
and only by rejecting these "false" frequencies 
could results of any reliability be obtained. The 
results are plotted on the usual diagram in figure 
17. Many of the points shown evidently are not 
related to the boundary-layer oscillations. 

2. Sound Through Hole in Plate 

A localized disturbance was produced by bring- 
ing sound into the boundary layer through a hole 
in the plate. At a position 18 inches above the 
center line and 4 feet from the leading edge a } -inch 
hole w T as drilled through the plate, and on the 
back was mounted a small headphone from a 
radio head set. This produced a three-dimensional 
disturbance emanating from a known location. 

A hot wire, mounted on the small sledlike head, 
was arranged for traversing upstream and down- 
stream at a fixed distance from the surface, usually 
along a line slightly displaced from the hole. A 
survey of possible types of measurement showed 
that the method could be used only for obtaining 
the type of result shown in figure 18. With the 
regeneration hook-up, that is, with the amplified 
output of the wire fed into the headphone, the 
frequencies of oscillation were determined with 
the wire at various distances downstream from 
the hole. The speed was constant during a run. 
Oscillation of the system occurred most easily for 
values of U less than 40 feet per second and with 
the hot wire less than 0.1 inch from the surface. 
The solid curves of figure 18 show the discon- 
tinuous manner in which the frequency varied 
with the distance. Reversing the connections to 
the headphone changed the phase of the input by 
180 degrees and gave the results indicated by the 
broken curves. 

The explanation of this sort of performance is 
that proper phase relations had to exist for re- 
generation and that the frequency automatically 
decreased as the distance was increased in order 
to keep an integral number of waves between the 
hole and the wire. There existed a preferred fre- 
quency band, determined by the characteristics 
of the wire, electrical and acoustic circuits, and 
boundary layer. Jumps occurred when another 



wave between hole and wire was added to keep 
the frequency as near as possible to the center of 
the preferred band. The distance between curves 
of one type along any line of constant frequency 
is thus the wavelength corresponding to that 
frequency. 

It is believed that the frequency band was- 
determined more by the boundary layer than by 
the rest of the circuit when the wire was far from 
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Points show regeneration frequencies resulting from feedback from hot wire 
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the hole and that the preferred frequency was 
that for which amplification was a maximum over 
the path from hole to wire. The frequency mid- 
way between the extremes and as far as possible 
from the hole was therefore read from sets of 
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curves for different wind velocities like that in 
figure 18 and placed on the usual Schlichting 
diagram. The result is shown in figure 19. The 
path over which amplification could take place is 
indicated by the horizontal line connected to each 
point. 

As the wavelengths corresponding to these fre- 
quencies were now known, it became possible for 
the first time to plot the results in terms of ad* 
and B. The theoretical curve and the experi- 
mental results are shown in figure 20. The lines 
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Figure 20. — Wavelength diagram showing zone of amplifi- 
cation according to Schlichting. 

Points show values of ad* for frequencies of figure 19. Lines from points show 
path of amplification. 

of total amplification for constant jft r „/£7 2 , corres- 
ponding to those in figure 19, are known. The 
paths over which amplification could occur in the 
experiment are represented as before by lines ex- 
tending from each point. Figures 19 and 20 
generally show the kind of results predicted by 
theory, if the frequency and wave length are de- 
termined mainly by the characteristics of the 
boundary layer. 

3. Excitation by Vibrating Ribbon 

(a) The Ribbon 

All attempts to excite boundary-layer oscilla- 
tions by wires vibrating in the boundary layer 
failed because the disturbance of the eddying 
Wake from the vibration of wires large enough to 
affect the layer produced early transition. Flat 
strips outside the layer suffered from self-induced 



flutter and also had troublesome wakes that, 
carried with them fluctuating pressure fields far 
downstream. A thin flat ribbon of phosphor 
bronze, placed edgewise to the flow in the bound- 
ary layer and only a few thousandths of an inch 
from the surface, was finally found capable of 
oscillating the boundary layer but produced no 
detectable effect when not vibrating. The ribbon 
was driven in and out from the surface by passing 
through it a small current of the desired frequency 
in the presence of a strong magnetic field from an 
electromagnet on the opposite side of the plate. 

Several ribbons were tried to find the best width, 
distance from surface, and length of vibrating 
span. All ribbons were 0.002 inch thick. Widths 
ranged from }{ 6 to }i inch; distances from surface,, 
from 0.006 to 0.015 inch; and vibrating segments, 
from 6 to 12 inches. None of these factors ap- 
peared to affect the performance to any noticeable 
extent. A certain amount of tension was neces- 
sary for stability. When it was desired to place 
the resonant frequency above the working range, 
high tension was necessary. 

After this preliminary work, a ribbon 0.1 inch 
wide and 3 feet long was installed 4 feet from the 
leading edge of the plate. At a distance of 6 
inches from each side of the center line of the plate, 
two strips of Scotch tape were laid on the surface 
for insulation and spacing. Another strip of tape 
was then laid on the ribbon to hold it firmly in 
contact with the tape on the plate, except for the 
segment 12 inches long at the center. Rubber 
bands at the ends supplied the tension. The 
12-inch segment at the center was 0.006 inch from 
the surface and was free to vibrate under the action 
of the current and the magnetic field. Examina- 
tions with a microscope showed that the ribbon 
vibrated to and from the surface in a single loop. 
The same oscillator and amplifier that had been 
used for the loudspeaker were found to furnish 
ample power for driving the ribbon. Different 
positions for the ribbon were chosen to cover a 
large range of Reynolds number, but the method 
of attachment was always the same. Distance 
from the surface varied somewhat for different, 
installations. 

(b) Procedure 

Two hot wires, one 0.010 inch and the other 
0.110 inch from the surface, both sensitive only to 
u, were mounted on the small sled-like head and 
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arranged for sliding upstream and downstream 
along the center line of the plate with the distance 
of the wires from the surface fixed. The fluctua- 
tions introduced into the boundary layer could be 
picked up by either hot wire. The wave form 
showed some distortion for an inch or so behind 
the ribbon but in general became nearly sinusoidal 
at greater distances. Frequencies amplified by 
the boundary layer showed increasing amplitude 
with increasing distance from the ribbon, while 
those damped showed decreasing amplitude. 
For much of the working region the oscillations 
could probably be regarded as two-dimensional, 
but no attempt was made to verify this assumption. 

Just as for the loudspeaker and the headphone 
behind the hole in the plate, the system could be 
made to oscillate by connecting the amplified out- 
put of the hot wire to the ribbon. Curves like 
those shown in figure 18 were obtained and, from 
them, three points plotted in each of figures 19 and 
20 were obtained. Since the ribbon could be used 
for applications that were believed to be more 
important, very little work was done with the 
regeneration hook-up. 

The wavelength was determined in the following 
simple manner: The input to the ribbon from the 
oscillator was connected to one pair of plates in the 
cathode-ray oscillograph and the output from the 
wire was connected to the other pair. A station- 
ary Lissajou figure consisting of a single closed 
loop was obtained since the frequencies of both 
input and output were the same. As the spacing 
between ribbon and hot wire was changed this 
figure changed from a straight line to an ellipse, 
then to a circle, again to an ellipse, and finally to 
a straight line inclined 90 degrees to the first line. 
This indicated that the phase between input to the 
ribbon and output from the wire had changed by 
180 degrees and that the change in spacing be- 
tween wire and ribbon was one-half wavelength. 
With the oscillograph serving as the indicator, 
wavelength measurements were thus reduced to a 
measurement of the distance moved by the hot 
wire. The w^ave velocity was then obtained by 
multiplying the wavelength by the frequency of 
the oscillations. 

Determinations of amplification and damping 
were carried out by setting the ribbon into 
vibration at an arbitrary fixed amplitude and 
reading the mean-square hot-wire output voltage 
as the wire was moved away from the ribbon. 



The results obtained are given in section VI, 3 (d). 
The distribution across the boundary layer was de- 
termined by traversing normal to the surface at 
some fixed distance behind the ribbon with a 
single hot wire sensitive to 11, mounted on the 
traversing apparatus described in section II, 7. 

During the course of the work frequencies 
ranging from 10 to 260 cycles per second were 
used. 

(c) Distribution of Amplitude 

Since a theoretical distribution of amplitude 
across the boundary layer has been given by 
Schlichting (see section V, 2 (d) and table 4), the 
actual distribution is of interest as a test of agree- 
ment between theory and experiment. Traverses 
across the boundary layer were therefore made 
with a hot wire downstream from the ribbon with 
frequency, air velocity, and distance from the 
leading edge of the plate chosen to conform to the 
conditions for which the theoretical calculations 
were made. These conditions correspond to 
oscillations on the neutral curve specified as 
follows: 

Conditions required by theory 
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1 Ribbon about 0.009 in. from surface. 

The distributions are given in figure 21 with 
amplitude expressed in terms of u'/U . The solid 
curves show the theoretical distribution and the 
broken curves through the points show the ob- 
served distribution. In both cases the traverses 
were made 4 inches behind the ribbon. The 
vibrational amplitude of the ribbon was set 
arbitrarily and about the same for branches I and 
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II. The smaller values of u'/U in branch II 
result from the smaller effect of the ribbon on the 
thicker boundary layer. Since the theoretical 
values given by Schlichting were based on an 
arbitrary average amplitude of u'/U = 0.05 from 
y=0 to y=8, Schlich ting's values were reduced 
to give the average found by experiment; that is, 
areas under the theoretical curves w^ere made 
equal to areas under tli3 experimental curves 
from y=0 to y=d. 

Although the general agreement between theory 
and experiment is good, perhaps the most striking 
confirmation of the theory is the 180 degree phase 
shift occurring at y 15 = 0.7. On moving outward 
from the surface all oscillations disappeared at 
this point and farther out they reappeared but 
with phase opposite to those that had disappeared. 
This phenomenon is shown in a more striking 
manner by the oscillograms in figure 22. The 
pairs of traces are simultaneous records of oscil- 
lations picked up by two wires, one of which was 




Figure 21. — Distribution of amplitude of oscillations 
across boundary layer. 

Solid curves are theoretical according to Schlichting. Broken curves are 
experimental. 
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Figure 22. — Oscillograms showing phase reversal in 
u-component of oscillations excited by vibrating ribbon. 

Simultaneous records from two hot wires located 1 ft downstream from 
ribbon. Lower trace— hot wire 0.055 in. from surface. Upper trace— hot 
wire at various distances from surface. Ribbon 3 ft from leading edge. 
Frequency = 70 cps. Uo=42 ft/sec. 

at a fixed distance of 0.055 inch from the surface 
(lower trace) and the other at varying distances, 
beginning at 0.080 inch (upper trace). The simul- 
taneous records were made possible by an electronic 
switch that impressed the separately amplified 
voltages of the two wires on the oscillograph, 
alternately and in rapid succession. The switching 
from one to the other accounts for the dotted 
appearance of the trace. The waves are in phase 
for the first four records; the amplitude of the 
upper wave decreases with each succeeding record 
to zero in the fifth. Beginning with the sixth, an 
oscillation reappears in the upper trace but 180 



Laminar Boundary-Layer Oscillations 



279 



degrees out of phase with that in the lower trace. 
The oscillations in these records lie near the 
center of the amplification zone at R= 1,720. 

The type of distribution shown in figure 21, 
including the phase reversal, was found for am- 
plified, damped, and neutral frequencies and for 
oscillations as close as 1 inch behind the ribbon 
and as far downstream as the oscillations remained 
free from transition disturbances. The fact that 
the main features of the distribution were found 
1 inch behind the ribbon shows that the forced 
oscillations in the boundary layer very soon reach 
a steady state. 

(d) Amplification and Damping 

When the wire was kept at a fixed distance from 
the surface and moved downstream from the 
ribbon, the amplitude of the forced oscillation 
either increased, decreased, or remained the same, 
depending on the frequency, the velocity, and the 
Reynolds number. The amount of the change 
with distance could be determined by reading the 
mean-square value of the fluctuating voltage 
across the wire on the output meter of the turbu- 
lence amplifier. In order to eliminate as far as 
possible the extraneous disturbances introduced 
by turbulence and by the noise of the tunnel, the 
signal from the wire was passed through an elec- 
trical filter, which was tuned in each case to pass 
only a 5-cycle band around the frequency being 
measured. The results of a typical run with the 
ribbon at the 4-foot position are shown in figure 
23. In this figure x denotes a position 2 inches 
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Figure 23. — Growth and decay of u-component of oscilla- 
tions produced by vibrating ribbon 4 ft- from leading edge 
of plate. 

xo is 2 in. downstream from ribbon. Uo=64 ft/sec. 



downstream from the ribbon. The ratio u'/u is 
the square root of the meter reading at x divided 
by the square root of the meter reading at x . 

According to eq 31 



or 



u' C t{x) 

—i=exv fitdt 
u r Jux,) t 

2.3 log 10 -r= 04t 

^0 Jt(xo) 



Differentiating with respect to t and using the 
relation dx/dt=c T> where c T is the wave velocity, 
gives 



0,=2.3c r 



( logi ° i) 



dx 



(39) 



The wave velocity is obtained by multiplying the 
measured wavelength by the frequency. If the 
changes in u'/u with distance shown in figure 23 
are truly effects of amplification and damping, 
the slopes of the curves in this figure serve to 
determine /?*. 

Before these slopes can be used with confidence, 
effects introduced by the thickening of the layer 
with increasing x and effects from possible changes 
in the distribution across the layer also with x 
must be investigated. Errors can arise from these 
sources because of the experimental procedure, 
namely, the determination of growth or decay of 
amplitude by varying the distance between the 
ribbon and the wire while keeping the wire at a 
fixed distance from the surface. From the evi- 
dence in section VI, 3 (c), changes in the distribu- 
tion are small and can cause but little error; how- 
ever, the increase of 5 with x may give rise to an 
apparent amplification or damping, depending on 
the position of the wire in the layer. In all cases 
in which data on amplification and damping were 
obtained, the wire was in the rising part of the 
distribution, that is, values of y/8 below those for 
maximum v'/U in figure 21. Consequently, the 
increasing 8 with increasing x introduced an ap- 
parent damping, which caused all observed ampli- 
fication to be less than the true value and all 
observed damping to be too great. If it is as- 
sumed that the inner part of the distribution curve 
is a straight line passing through the origin, the 
magnitude of the error in f} t calculated from eq 
39 by use of observed slopes is (%)c T /x. Any cur- 
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vature in the distribution curve like that shown 
between zero and the peak value of u'/U in figure 
21 lessens the error. Since (%)c r /x therefore 
represents an upper limit, and since its value 
turned out to be about the order of the random 
deviation in results, p t was calculated by eq 39 
with the use of the experimentally determined 
slopes without correction. 



and the corresponding values of I\. This amounts 
to finding points of zero slope on curves like those 
shown in figure 23. Such points were obtained 
from a number of runs at various wind speeds 
with the ribbon at the 4-foot position and again 
with the ribbon 9 inches from the leading edge. 
Values of p r v/U 2 and the corresponding values of 
R were computed and plotted in figure 24. The 
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Figure 24. — Frequency of neutral oscillations excited in boundary layer by vibrating 

ribbon. 

Solid curve is neutral curve according to Schlichting. Broken curves are neutral curves denned by experimental 

points. 



In figure 23 the curve for 40 cycles per second 
shows damping throughout. As the frequency is 
increased amplification takes place, reaches a max- 
imum at 120 cycles per second, and then decreases 
until 180 cycles per second again shows damping. 
In the range of Reynolds number represented by 
this range in x, all frequencies below about 50 and 
above 180 cycles per second are damped, while a 
band of amplified frequencies lies in between. 

Of considerable interest are the frequencies for 
which there is neither amplification noi* damping 



solid curve is again Schlich ting's neutral curve 
and the broken curves are those defined by the 
experimentally determined points. The agree- 
ment is very satisfactory, except for points above 
^/c7 2 =200X10- 6 . 

By use of the measured wavelengths, ad* was 
calculated for the neutral frequencies and the 
results were plotted on the wavelength diagram 
shown in figure 25. The open circles denote 
points on branch II, and the closed circles denote 
points on branch I. The broken lines drawn 
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through these points define the experimental 
neutral curves. Schlichting's theoretical curve 
is shown and also the theoretical neutral curve 
given by C. C. Lin. 4 

Still another type of representation involving 
both frequency and wavelength is the diagram of 
c r /U and R in figure 26. The points and the 
broken curves show the wave velocities found for 
the neutral oscillations. The solid curve is 
Schlichting's theoretical curve. 

It is again pointed out that both branches of 
the theoretical curves in figures 24 to 26 meet the 
i?-axis at R= oo and thus extend beyond the 
region indicated. Transition prevented experi- 
mental determinations above R= 2,800. 

.5 



The wavelength diagram is repeated in figure 27 
and this time shows, in addition to the theoretical 
neutral curve, the curves of equal amplification 
within the zone. These contours were copied 
directly from a figure given by Schlichting [4]. 
The broken curves are defined by the neutral 
points, which are not shown in this figure. The 
points with the values opposite (all to be multi- 
plied by 10~ 4 ) are the observed values of pi5*/U 
obtained from values of the slopes other than 
zero. Negative values signify damping. 

The distribution across the amplification zone 
at particular values of R is shown in figure 28. 
The theoretical curves, drawn solid, are plotted 
from values read from Schlichting's figure in [4]. 
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Figure 25. — Wavelength of neutral oscillations excited in boundary layer by vibrating 

ribbon shown by broken curves and points. 

Closed circles— branch I. Open circles— branch II. Theoretical neutral curves shown by solid lines. 

The broken curves are drawn through the experi- 
mental points. 

The comparison between theory and experi- 
ment given in figures 24 to 28 shows, on the whole, 
better agreement than might be expected when 
it is considered that at the outset the very exis- 
tence of such phenomena had not been proved. 
The observed amplification shown in figures 27 
and 28 is less than the theoretical, but it will be 
recalled that the theoretical calculation yielded 
accurate values only near the neutral curve. In 
the neighborhood of the low Reynolds number 
limit the experimental inaccuracies are large be- 
cause of the small amplification and damping. 
The maximum values of /3 r v/U 2 , a<5*, and c r /Tj 
and the minimum value of R for which amplifica- 
tion can occur as determined approximately by 
experiment are compared with the theoretical 
limits in the following table: 
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Figure, 26. — Wave velocity o] neutral oscillations excited 

in boundary layer by vibrating ribbon. 
Solid curve is neutral curve according to Schlichting. Broken curves are 

neutral curves defined by experimental points. Closed circles— branch I. 

Open circles — branch II. 

4 Lin's curve was not available when this work was originally published. 
The curve is a recent addition taken from C. C. Lin, On the stability of two 
dimensional parallel flows, Quarterly of Applied Mathematics, III, [2], 117 
to 142 (July 1945); [3], 218 to 234 (Oct. 1945); [4], 277 to 301, (Jan. 1946). 
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Tollmien's limiting values, on the whole, agree 
better with experiment than Schlichting's values. 
A comparions of the experimental results in figure 
25 with Tollmein's diagram of a8* plotted against 
log 10 R in [2] shows better agreement in general 
than with the corresponding diagram given by 
Schlichting in figure 25. In spite of this, Schlicht- 
ing's results were used because they were more 
complete than Tollmien's. 
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Fic.tjre 27. — Contours of equal amplification according to Schlichting. 
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Figure 28. — Distribution of amplification across zones of 

figure 27. 

Solid curves are theoretical. Broken curves are experimental. 

VII. Effect of Pressure Gradient on 
Boundary-Layer Oscillations 

Figure 29, obtained in connection with the 
earlier work presented in section V, 1, shows the 
effect of pressure gradient on naturally excited 
oscillations. The pressure distribution along a 
part of the plate is shown to the left in this figure. 
The reduction of amplitude in the region of falling 
pressure is very marked. Accelerated growth of 
amplitude is evident in the region of rising pres- 
sure. It was accordingly inferred that negative 
pressure gradients decreased amplification (or 
increased damping) of the oscillations while posi- 
tive gradients increased amplification. A quan- 
titative investigation of this effect was therefore 
undertaken with oscillations forced on the bound- 
ary layer by the vibrating ribbon. 

When the ribbon was in the 3-foot position, the 
regular procedure of measuring amplifications 
and damping and of mapping the neutral curves 
was carried out for pressure distributions B, C, 
D, and E shown in figure 30. The corresponding 
pressure gradients are given in table 5. 

In figure 30, q x is the dynamic pressure at any 
distance x from the leading edge, and q is a con- 
venient reference pressure. Curve A shows one 
of the distributions for the condition generally 
termed "zero pressure gradient." This condition 
was obtained by giving sufficient divergence to 
the adjustable side walls of the working chamber 
to overcome the normal pressure drop. The 
shape of the curve in the first 6 inches was deter- 
mined by the flow around the leading edge of the 



plate and hence could be changed but little by 
the setting of the walls. The walls also had 
little effect on the small local variations, since 
these were found to be produced by slight wavi- 
ness of the surface. Curve A is thus a condition 
of zero gradient only when the length of the 
interval and the end points are properly chosen. 
With a pressure distribution given by curve A, 
the experimental values of velocity distribution 
and 5* were in good agreement with the Blasius 
values; hence, the approximation to zero gradient 
is probably sufficiently good. 



-.05 -.10 











Figure 29. — Oscillograms of u- fluctuations showing effect 

of pressure gradient on boundary layer oscillations. 

Distance from surface=0.021 in. L T o=95 ft/sec. Time interval between 

dots =1/30 sec. 

Curve B, which shows the pressure rising along 
the plate, was obtained by giving the walls their 
greatest possible divergence. Curve C, which 
shows the pressure falling, was obtained with the 
walls forming a converging passage. The steep 
positive gradient shown by curve D was produced 
by a vertical airfoil placed 10 inches from the side 
of the plate running from the floor to the ceiling 
of the tunnel, and the steep negative gradient of 
curve E was produced by blocking the stream 
with a plywood panel % inch thick and 28 inches 
wide extending from the floor to the ceiling placed 
parallel to and 10 inches from the plate. 
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The results obtained on the effect of pressure 
gradient are given in figures 31 to 34 for distribu- 
tions B and C and will be discussed in detail later 
in this section. Corresponding figures for dis- 
tributions D and E could not be obtained because 
of the large effect of these gradients on the oscil- 
lations. With the large gradients, the general 
impression given by observing the oscillations was 
that amplification always occurred in gradient 
D regardless of the frequency and Reynolds num- 
ber, while damping always occurred in gradient 
E\ however, this statement must be somewhat 
qualified. For example, in E fluctuations of all 
frequencies were damped for values of R up to 
2600, which w T as the highest value obtainable with 
the ribbon at the 3-foot position. It is quite pos- 
sible that a zone of amplification lay beyond the 
range of the experiment. In the case of gradient 
D, the investigation had to be restricted to low 
Reynolds numbers because of the occurrence of 
transition. The few observations made here indi- 
cated a possible low limit of Reynolds number 
below which damping occurred, but the results 



are in doubt because of a possible scale effect on 
the airfoil used to produce the pressure gradient. 
The pressure rise may have been lessened by lami- 
nar separation on the airfoil at the low velocities 
and the positive gradient may have become less 
than that given by D. Accurate measurements 
of pressure distribution at the low velocities was 
difficult. 

Before presenting the results for the two small 
gradients B and C, attention must be called to the 
effect of pressure gradient on <5*. Since 5* enters 
into certain dimensionless quantities used in plot- 
ting the results, its actual value must b.e known for 
each case. The value of 5* was measured only for 
distribution B and the result was 



Observed 5* 

Hhisius ,5* 



= 1.09 at z=3.33 ft. 



For distribution C, <5* was estimated by the aid of 
the Karman-Pohlhausen method by assuming a 
linear pressure distribution. The validity of the 
met hod was tested by making a similar calculation 
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Figure 31. — Frequency diagram of neutral oscillations. 
Experimental curve for zero pressure gradient shown without points. Posi 
tive gradient from distribution B and negative gradient from distribution 
C of table 5. 
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Figure 32. — Wavelength diagram of neutral oscillations. 

Experimental curve for zero pressure gradient shown without points. Posi- 
tive gradient from distribution B and negative gradient from distribution 
C of table 5. 
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Figure 33.— Wave-velocity diagram of neutral oscillations 
Experimental curve for zero pressure gradient shown without points. Posi- 
tive gradient from distribution B and negative gradient from distribution 
C of table 5. 
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Figure 34. — Distribution of amplification across zones of 
figure 31, showing effect oj pressure gradient. 

Positive gradient from pressure distribution B and negative gradient from 
distribution C of table 5. 

for B and obtaining good agreement with the 
measured 5*. The estimated result for C was 



Calculated 3* 
Blasius 5* 



=0.93 at z=3.33 ft. 



The values of 5* for the Blasius distribution were 
multiplied by the above factors to give the value 
of 5* pertaining to gradients B and C. 

Figure 31 is the much-used frequency diagram 
showing the results for zero, positive, and negative 
gradients of distributions A, B ) and C, respectively. 
The solid curve is the same as the broken curve of 
figure 24 for zero gradient, shown here without the 
experimental points. The points for negative 
gradient lie close to this curve; the points for 
positive gradient generally lie outside the zone of 
amplification for zero gradient. The positive 
gradient has therefore widened the amplification 
zone as expected. Contrary to expectation, the 
negative gradient failed to reduce the size of the 
zone. No observations were made below i? = 800. 

Since wavelengths of the oscillations were 
measured, the data for these same frequencies are 
given in a wavelength diagram (fig. 32) and a 
wave-velocity diagram (fig. 33). It will be noted 
that figure 32 is unlike the usual wavelength 
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diagram (for example, fig. 25) in that the ordinate 
is now av/U instead of a<5*. This change was 
made in order to eliminate 5* from one of the 
coordinates and thereby to restrict errors occurring 
in the estimated value of 5* with pressure gradient 
to the Reynolds number coordinate. The new 
ordinate is simply a<5* divided by the Reynolds 
number, that is, 

a5* __ av 



In both figures 32 and 33, the solid curves are the 
experimentally determined neutral curves for zero 
gradient with the experimental points omitted. 

Figures 32 and 33, like figure 31 , show an effect 
of the positive gradient but little or no effect of 
the negative gradient. These figures show only 
the neutral points and so merely define the bound- 
aries of the amplification zone. The effect of 
pressure gradient on the amplification within the 
zone is demonstrated in figure 34. This figure, 
which is similar to figure 28, except that zones are 
from the frequency diagram instead of the 
wavelength diagram, shows the distribution of 
amplification coefficient, actually Pi8*/U , across 
the zones of figure 31 at certain fixed values of the 
Reynolds number. The positive and negative 
pressure gradients signify as before pressure dis- 
tributions B and C, respectively. For compari- 
son, the experimental curves for zero gradient arc 
given as solid lines. The comparision is best made 
between pairs of curves having Reynolds numbers 
most nearly alike. 

Since P^/Uo increases with Reynolds number 
(see fig. 27), the difference between R= 1,840 for 
zero gradient and if =1,810 for negative gradient 
may account for the reduced amplification with 
negative gradient. The failure to show conclusively 
an effect of the small negative gradient may be 
due in part to the small local value of the gradient 
shown in table 5. The difference between R= 
2,200 for zero gradient and R= 2,050 for positive 
gradient is significant, since R= 2,050 for zero 
gradient would lie below R= 2,200. 

VIII. Transition From Oscillations to Tur- 
bulence, Zero Pressure Gradient 

The distinguishing feature of transition with low 
stream turbulence is the occurrence of boundary- 
layer oscillations that grow until the layer becomes 



turbulent. Under such conditions, transition to 
turbulent flow is simply transition from a particu- 
lar type of oscillatory motion to turbulent motion. 
The oscillations are now understood and (he final 
step in an understanding of the origin of turbulence 
involves an investigation of this transition process. 

Table 5. — Pressure gradientb (l/qi) (AP/A.r) 1 

(/i=dynamic pressure outisde boundary layer at 3.33 n from leading edge of 

plate. 
P=pressure near surface at distance i from leading edge. 



Distribution 


Average 

gradient 

from 1 to 

6 ft 


Local gradient 


A 

B 

C. . 


Per edit 
-0. 11 
1 . 62 
-2. 38 


Percent 
-0.78 from 4 to 5 ft 

2. 48 from 3 to 3.67 ft 
—1.36 from 3 to 3.67 It 


D 


12. from 3 to 3.67 ft 


E 




—9. 8 from 3 to 3.67 ft 









1 Note that pressure gradient and dynamic pressure gradient have opposite 
signs. (See fig. 30.) 

Although transition in highly disturbed flow 
does not fall into this category, the two are similar 
in that both have to do with the effect of a dis- 
turbance large enough to bring about turbulent 
flow. The effect, however, depends on the type 
of the disturbance as well as its magnitude; and, 
as has been shown, an amplified disturbance does 
not have the random character of an initially 
large disturbance. Transition from an amplified 
disturbance promises to be of considerable im- 
portance in modern low-turbulence wind tunnels 
and in the atmosphere. 

The vibrating ribbon to produce oscillations 
and the hot-wire apparatus to examine their 
development were suitable tools for the investi- 
gation of transition. The ribbon was placed at 
the 3-foot position and two hot wires at different 
distances from the surface were arranged for tra- 
versing parallel to the plate through the transition 
region. The Reynolds number was always chosen 
so that the undisturbed boundary layer (ribbon 
not vibrating) was laminar at the hot wires. 
Representative samples of the results obtained 
are shown in figures 35 to 37. 

Figure 35 shows what happens at a fixed point 
2.5 feet downstream from the ribbon as the vi- 
brational amplitude of the ribbon is varied. The 
pairs of traces show simultaneously the response 
of the two wires. The values of gain control 
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Fi-gure 35. — Oscillograms of u-fluctuations showing devel- 
opment of turbulence from oscillations produced by vibrat- 
ing ribbon at position 2.5 ft down-stream from ribbon, 
Kibbon 3 ft from leading edge of plate. Simultaneous records at 0.014 in- 
born surface, upper trace; and 0.114 in. from surface, lower trace. Fre" 
quency 70 cps. U =50 ft/sec, 

setting on driving amplifier are roughly propor- 
tional to ribbon amplitude. The amplitude of the 
waves on the traces cannot be quantitatively re- 
lated to fluctuations of velocity because the sensi- 
tivity of the wires was unknown. The sensitivi- 
ties of the two wires were different, and so com- 
parison of amplitudes at inner and outer wires was 
impossible. It is known, however, that velocity 
increases produce upward displacement of the 
traces and that velocity decreases produce down- 
ward displacement. The direction of time in this 
figure is from left to right. 

This figure shows that the boundary layer at a 
given point can be changed from laminar to turbu- 



lent by varying ribbon amplitude and shows 
qualitatively the progressive nature of the change. 
Examination of successive pairs of traces indicates 
that fairly regular waves first make their appear- 
ance, being first in phase at the two wires and 
later in the fifth and sixth records, in and out of 
phase in random fashion as the amplitude in- 
creases. The distances from the surface are such 
that an in-phase condition would be expected 
from the distribution of figure 21. A growing 
distortion of the waves then appears until, in 
the seventh record, bursts of high-frequency fluc- 
tuations occur. The latter, identified by the 
jagged appearance of the trace, are bursts of 
turbulence. The final record shows the com- 
pletely turbulent state. The dotted character of 
the. traces, due to the electronic switch, obliter- 
ates details of rapid fluctuations. A similar suc- 
cession of changes is shown in figures 36 and 37. 
In figure 36 simultaneous records are again 
shown, this time from two wires on the same 
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Figure 36. — Oscillograms of u-fluctuations showing devel- 
opment of turbulence from oscillations produced by vibrat- 
ing ribbon placed 3 ft from leading edge of plate. 

Simultaneous records at 0.011 in. from surface, upper trace; and 0.112 in. from 
surface, lower trace. Frequency 70 cps. U Q =50 ft/sec. 



288 



Journal of Research 



mounting moved downstream from the ribbon 
while the amplitude of the ribbon remained 
fixed. Figure 37 shows the same condition but 
for the two traces photographed separately in 
order to eliminate the dotting effect of the elec- 
tronic switch. The difference in the appearance 
of the records in these two figures is due more to 
the different magnification than to the character 
of the traces. Since the frequency of 70 cycles 
per second is amplified by the boundary layer, 

.Oil IN. FROM SURFACE 



velocity at the hot wire. The percentages are the 
same for all records in figure 36. 

For fluctuations as large as most of those 
shown, it is not possible to indicate the zero- or 
mean-velocity position of the trace. This is 
because of the nonlinear characteristics of the hot 
wire which produces greater displacements of the 
trace for velocity decreases than for increases. A 
certain amount of distortion and some error in 
the records will result from this effect. It is 
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Figure 37. — Oscillograms of u- fluctuations shoiving development of turbulence from oscillations produced by 
vibrating ribbon placed 3 ft from leading edge of plate. 
Frequency 70 cps. U =50 ft/sec. 



the amplitude increases downstream and increas- 
ing distance is in one respect like increasing ribbon 
amplitude in figure 35. The two cases are not 
equivalent, however, because the boundary-layer 
thickness increases from record to record in 
figures 36 and 37. 

From calibration data, obtained in connection 
with figures 36 and 37, it was possible to relate 
two particular displacements on the record to a 
percentage change in the mean velocity at the 
position of the wire. These two percentages are 
indicated at the right of each trace, an upward 
displacement indicating an increase and a down- 
ward displacement a decrease from the mean local 



certain, however, that crests are above the mean 
velocity while troughs are below the mean and 
the distance between the two is approximately the 
double amplitude. The scale at the right fur- 
nishes a measure of the double amplitude. 

In the succession of changes occurring from 
point to point through the transition region, the 
oscillations first increase in amplitude and, once 
large, undergo increasing amounts of distortion. 
Turbulence appears at first in short intervals of 
the trace and finally for the full length of the trace. 
Randomness occurs gradually, and this suggests 
a continuous transition from wave to turbulent 
motion. It is not difficult to imagine a process 
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here like that often assumed for the formation of 
eddies from a free vortex sheet. The sheet is 
imagined to take first a wavelike character, then 
as the wave grows, to curl up into discrete eddies. 
The disturbed laminar boundary layer may be 
regarded as a wavy vortex layer with the wave 
progressively increasing in amplitude and dis- 
torting until discrete eddies are formed. The 
eddies themselves are unstable and soon break up 
into a diffusive type of motion which characterizes 
turbulent flow. If the validity of this picture is 
granted, transition is understood when the wave 
motion is understood. An understanding of 
transition is therefore close, but only close because 
the wave motion is understood only when the 
amplitude is small. Furthermore, the important 
question, concerning where along the surface 
eddies first form, is still unanswered. 

Another process may also be imagined. It will 
be noted in figures 35 to 37 that transition takes 
place in a more orderly manner at the greater 
distance from the surface. Here turbulence 
appears in certain parts of the wave cycle, gener- 
ally in the sharp-pointed low-velocity part of the 
cycle. This behavior lends support to the view 
held by some investigators that transition occurs 
when the flow comes to rest or reverses direction 
near the surface (separation point). While the 
amplitudes shown do not reach the necessary 100 
percent, it seems possible that this value would 
be found at some point nearer the surface. Against 
this argument is the fact that the records made 
nearer the surface do not show the turbulent 
bursts confined to the low-velocity side of the 
cycle. The argument is further weakened by the 
possibility of accounting for the observed effects 
at the greater ?/- values by the difference between 
the distribution of mean-velocity characteristic of 
a laminar layer and that characteristic of a 
turbulent layer. This is to say, that the observed 
effect is caused by the turbulence rather than a 
cause of the turbulence. 

The idea of transition by separation has been 
formulated by Taylor [1] by assuming inter- 
mittent separation or near separation to be caused 
by the local pressure gradients accompanying 
disturbances. The assumption is that separation 
is imminent when the Karrran-Pohlhausen para- 
meter A, derived from pressure gradients either 
steady or fluctuating, reaches a certain negative 
value. A steady adverse pressure gradient large 



enough to cause laminar separation is known to 
bring about transition rather than complete 
separation when the Reynolds number is suffi- 
ciently high. By analogy, intermittent separa- 
tion is supposed to have the same effect. For 
isotropic turbulence in the surrounding flow the 
fluctuating pressure gradient is known in terms of 
measurable quantities, namely, u' and the scale. 
Taylor's relation is 






i(iT -*<*>> 



where L is the scale, and F (R x ) expresses a func- 
tion of the ^'-Reynolds number of transition. 

In essence Taylor's relation rests on the concept 
of velocity fluctuations impressed on the boundary 
layer from pressures outside. When velocity 
fluctuations belong to the natural wave in the 
boundary layer, Taylor's relation can have little 
bearing on the problem. There is therefore no 
theory connecting oscillations to transition. The 
idea of transition by separation leads to the con- 
cept that transition results when the 'U-component 
of the wave gets sufficiently large to reverse the 
flow near the surface. A lack of knowledge re- 
garding large-amplitude waves results in a lack of 
information concerning the point where transition 
will occur. 

It has not been proved that separation is a neces- 
sary condition for transition. Figure 11, where a 
slight reversal near the surface is indicated by a 
frequency doubling in the lower part of the cycle 
in the trace at 4 feet, gives evidence that reversal 
alone is not a sufficient condition. Perhaps the 
former picture, based on the analogy to the vortex 
sheet is the more accurate one. 

Since the Taylor relation is not applicable when 
transition results from amplified waves, the 
average scale of the initial disturbance is impor- 
tant only insofar as it is related to the spectrum of 
the disturbance. The frequency distribution of 
the initial disturbance is of great importance, 
since amplification finally selects only a very 
narrow band of frequencies from all those initially 
present to bring about transition. The presence 
or absence of frequencies near this amplified band 
is obviously more important than the total energy 
of the disturbance. 

Only on the basis of certain prominent frequency 
bands in the stream turbulence can figure 7 be ex- 
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plained. The reasoning is as follows: If no disturb- 
ances were present, theoretically transition should 
not occur, The pairs of curves in figure 7 should 
then approach R x = °o for zero turbulence. On 
the basis of amplification of initial distrubances in 
which all frequencies are present, equation 33 
suggests that the curves defining the transition 
region should fall rapidly at first and then at a 
progressively diminishing rate as the turbulence 
is further increased. Actually, the opposite trend 
is shown in figure 7. Figure 8 shows that not all 
frequencies are represented alike and, furthermore, 
that noise is responsible for many of the peaks 
shown. As the present turbulence was decreased 
by the addition of damping screens, the ratio of 
noise to true turbulence was increased. While 
the total energy of the distrubance was decreased, 
the energy in certain frequency hands may ac- 
tually have been increased. A reduction of 
noise might have greatly altered figure 7. 

IX. Conclusion 

Laminar boundary-layer oscillations are the 
velocity fluctuations that result from a wave 
traveling downstream through the boundary layer. 
The characteristics of this wave have 'been meas- 
ured and good agreement with theory has been 
obtained. The wave may be set up either by a 
vibrating object in the boundary layer, such as a 
vibrating ribbon, or by disturbances from the 
outside, such as stream turbulence and sound. 
When the disturbance is random or nearly ran- 
dom and contains many frequencies, amplification 
and damping isolate a wave containing a narrow 
band of frequencies in the neighborhood of the 
frequency most highly amplified. 

A wave in the boundary layer constitutes a 
disturbance that will cause transition on a flat 
plate with zero pressure gradient when the ampli- 
tude is sufficiently large. A sufficiently large 
random disturbance will also cause transition. 
When random initial distrubances are so small 
that transition is delayed until values of R^ 2,000 
have been reached, sufficient amplification can 
occur from the experimentally established lower 
limit R=450, to R^ 2,000 to give rise to a recog- 
nizable oscillation. Larger initial disturbances 
obscure the oscillations by irregular fluctuations 
of comparable magnitude. 

Natural disturbances may be sound waves as 
well as turbulence. Turbulence generally has a 



more nearly random distribution of energy with 
frequency than sound. In view of the importance 
of frequency, sound disturbances with a concen- 
tration of energy in frequency bands that are 
highly amplified by the boundary layer may he 
more conducive to early transition than turbu- 
lence. This is important in free flight where 
turbulence is probably negligible but where engine 
and propeller noise is present in large 4 amounts. 

The effect of pressure gradient on the oscillations 
is a practicably important phase of the problem 
about which more should be known. Little more 
has been done here 4 than to show the 4 direction of 
the effect and to give some idea of its magnitude. 
The results are in accord with theory, which also 
is incomplete in this respect. 

Other aspects of the problem remain to be in- 
vestigated. One of these having practical impor- 
tance is the effect of curvature. The importance 
of boundary-layer oscillations on airfoils cannot 
be determined until effects of curvature as well 
as pressure gradient are known. Other quantities 
associated with the wave hut not investigated are 
r and uv. A study of these quantities would be 
interesting as a further check on the theory, 
although hardly necessary to identify the oscilla- 
tions and their characteristics in view of the infor- 
mation derived from studies of u. 

It is possible that boundary-layer oscillations 
may arise from internal disturbances as well as 
from external disturbances — that is, from surface 
irregularities and vibration of the surface. A 
randomly distributed small roughness may pro- 
duce effects similar to small amounts of turbulence 
in the air stream. Small ridges or w r aves in the 
surface may start oscillations when the spacing is 
near some amplified wavelength. Vibration of 
the surface, like sound, may produce oscillations, es- 
pecially when the frequency is near some amplified 
oscillation frequency. An investigation of these 
and other phases of the problem will throw addi- 
tional light on the important problem of transition. 
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